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We present the BOOMERanG -03 experiment and maps of the Stokes parame- 
ters I, Q, U of the microwave sky obtained during a 14 day balloon flight in 2003. 
Three regions of the southern sky were surveyed: a deep survey (~ 90 square 
degrees) and a shallow survey (~ 750 square degrees) at high Galactic latitudes 
(both centered at RA ~ 5.5h, dec ^ —45° ) and a survey of ~ 300 square de- 
grees across the Galactic plane at RA ~ 9 Ah, dec ~ —47". All three surveys 
were carried out in three wide frequency bands centered at 145, 245 and 345 
GHz, with an angular resolution of ~ 10'. The 145 GHz maps of Stokes I are 
dominated by Cosmic Microwave Background (CMB) temperature anisotropy, 
which is mapped with high signal to noise ratio. The measured anisotropy pat- 
tern is consistent with the pattern measured in the same region by BOOMERanG 
-98 and by WMAP . The 145 GHz maps of Stokes Q and U provide a robust 
statistical detection of polarization of the CMB when subjected to a power spec- 
trum analysis. This amplitude of the polarization is consistent with that of the 
CMB in the ACDM cosmological scenario. At 145 GHz, in the CMB surveys, 
the intensity and polarization of the astrophysical foregrounds are found to be 
neghgible with respect to the cosmological signal. At 245 and 345 GHz we detect 
ISD emission correlated to the 3000 GHz IRAS/DIRBE maps, and give upper 
limits for any other non-CMB component. We also present intensity maps of 
the surveyed section of the Galactic plane. These are compared to monitors of 
different interstellar components, showing that a variety of emission mechanisms 
is present in that region. 

Subject headings: cosmology, cosmic microwave background, polarization, inter- 
stellar matter, polarimeters 

1. Introduction 

The Cosmic Microwave Background (CMB) is a remnant of the early Universe. Its 
existence is one of the pillars of the current Hot Big Bang model; its spectrum, temper- 
ature anisotropy, and polarization carry information about the fundamental properties of 
the Universe. The power spectrum of the temperature anisotropy of the CMB, (TT), is 
characterized by a flat plateau at scales larger than the horizon at recombination {9 » 1°; 
£ « 200), where primordial perturbations froze early in the history of the Universe, and 
by a series of peaks and dips at sub-horizon scales: the signatures of acoustic oscillations 
of the primeval plasma. Measurements of angular power spectrum have been very effective 
in constraining cosmological parameters (e.g., de Bernardis et al. 1994; Bond et al. 1998, 
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2000; Dodelson and Knox 2000; Tegmark and Zaldarriaga 2000a,b; Bridle et al. 2001; 
Douspis et al. 2001; Lange et al. 2001; Jaffe et al. 2001; Lewis and Bridle 2004; Netter- 
field et al. 2002; Ruhl et al. 2003; Spergel et al. 2003; Bennett et al. 2003; Tegmark et 
al. 2004). However, the temperature anisotropy power spectrum is degenerate in some of 
these parameters; independent cosmological information is required to break the degeneracy 
(Efstathiou and Bond 1999). All these studies make specific assumptions on the type of 
the initial conditions (adiabatic, or isocurvature) and on the shape of the power spectrum 
of the initial perturbations (power-law, scale- invariance, running index, etc.). When such 
assumptions are relaxed, the determination of the cosmological parameters becomes much 
more uncertain (e.g., Bucher et al. 2002). 

There is additional information encoded in the linear polarization properties of the CMB. 
CMB photons are last scattered at z ~ 1100. In Thomson scattering, any local quadrupole 
anisotropy in the incoming photons creates a degree of linear polarization in the scattered 
photons. The main term of the local anisotropy due to density (scalar) fluctuations is dipole, 
while the quadrupole term is much smaller. For this reason the expected polarization is quite 
weak (Rees 1968; Kaiser 1983; Hu & White 1997; Kamionkowski et al. 1997; Zaldarriaga 
2003). The polarization field can be expanded into a curl- free component (E-modes) and a 
curl component (B- modes). Six auto and cross power spectra can be obtained from these 
components: (TT), (TE), {EE), (BB), (TB), and (EB). Due to the parity properties of 
these components, standard cosmological models have {TB) — and (EB) = 0. Linear 
scalar (density) perturbations can only produce E-modes of polarization (e.g., Seljak et al. 
1997). In the concordance model, (EE) ~ 0.01 making (EE) a very difficult observable 
to measure. Tensor perturbations (gravitational waves) produce both E-modes and B-modes. 
If inflation happened (see, e.g., Mukhanov & Chibisov 1981; Guth & Pi 1982; Linde 1983; 
Kolb & Turner 1990), it produced a weak background of gravitational waves. The resulting 
level of the B-modcs depends on the energy scale of inflation, but is in general very weak 
(see, e.g., Copeland et al. 1993; Turner 1993). Alternative scenarios, like the cyclic model 
of Steinhardt & Turok (2002), do not produce B-modes at all (Boyle et al. 2004). 

Sensitive measurements of the polarization spectra will provide a conflrmation of the 
current scenario of acoustic oscillations in the early universe and improve the determination 
of cosmological parameters, in particular those related to the optical depth and reionization 
(see e.g. Kaplinghat et al. (2003)). They will help also in detecting deviations from a simple 
power-law spectrum of the initial perturbations. Moreover, they will allow study of the 
detailed mix of adiabatic and isocurvature initial perturbations (e.g., Gordon Sz Lewis 2003; 
Peiris et al. 2003). The detection of (BB) will probe the gravitational lensing of E-modes 
(Zaldarriaga & Seljak 1998), and, if present, the inflation generated component (e.g. Leach 
& Liddle 2003; Song & Knox 2003). 
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Confusion by galactic foregrounds will ultimately limit the precision with which {EE) 
and {BB) can measured. Not much is known about the galactic polarized background at 
microwave frequencies. The two mechanisms producing diffuse brightness of the interstellar 
medium are synchrotron radiation from relativistic electrons and thermal emission from dust. 
The former is sampled at low frequencies. Patchy high latitude observations at frequencies 
between 0.408 and 1.411 GHz are collected in Brouw & Spoelstra (1976). New observations 
carried out with the ATCA telescope at 1.4GHz (Bernardi et al. 2003) and at 2.3 GHz 
(Carretti et al. 2005) in the same high galactic latitude region observed by our experiment 
show that the polarized synchrotron emission is very weak. A naive extrapolation to 145 GHz 
predicts a polarization of 0.2 jiK rms, small with respect to the {EE) ~ AjiK expected 
in the concordance model. Polarized emission of galactic dust has been detected in the 353 
GHz survey of Archeops in the Galactic Plane (Benoit et al. 2003) and at high galactic 
latitudes (Pointhieu et al. 2005). There, {TE) has been detected at a level of 2 o", while only 
an upper limit was obtained for {EE). The polarized dust emission yH^EE) extrapolated 
to 145 GHz is quite weak, less than 1 /xK rms. While the foreground signals are expected 
to be smaller than the CMB signal at 145GHz, they should not be ignored for future, very 
precise measurements of CMB polarization. To do this, multiband measurements will be 
mandatory. 

After a long pioneering phase (Caderni et al. 1978; Nanos 1979; Lubin & Smoot 1981; 
Masi 1984; Partridge et al. 1988; Netterfield et al. 1995; WoUak et al. 1997), the measure- 
ment of CMB polarization is today a rapidly growing field; new interest has been sparked 
especially by the possibihty of detecting the {BB) signature of the inflationary gravity wave 
background (Keating et al. 2001; Subrahmanyan et al. 2002; Hedman et al. 2002; Piccirillo 
et al. 2002; Delabrouille et al. 2002; Masi et al. 2002; Villa et al. 2002; Kovac et al. 2002; 
Johnson et al. 2003; Keating et al. 2003; Kogut et al. 2003; Farese et al. 2004; Leitch et al. 
2004; Barkats et al. 2004; Rcadhcad ct al. 2004; Cortiglioni ct al. 2004; Cartwright ct al. 
2005). To date, statistically significant detections of CMB polarization have been reported 
by the experiments DASI, CAPMAP, CBI and WMAP, all using coherent techniques. DASI 
has detected {TE) at 2.9 a and {EE) at 6.3 a (Leitch et al. 2004); CAPMAP (Barkats et al. 
2004) has detected {EE) at ~ 2a; CBI (Readhead et al. 2004) has detected {EE) at > 7a; 
WMAP has detected {TE) at many a (Kogut et al. 2003). The polarization power spectra 
measured by these experiments are all consistent with the forecast from the "concordance" 
model best fitting the WMAP (TT) power spectrum. Their precision, however, is not yet 
good enough to improve significantly the constraints on the cosmological parameters. The 
only exception is the {TE) measurement by WMAP at large angular scales, which provides 
evidence for an early, complex reionization of the universe (Kogut et al. 2003; ?). The 
detailed structures in the {EE) spectrum are still to be confirmed, and we are very far from 
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the sensitivity required to constrain the initial conditions or inflation. 

The CMB polarization signals are so small with respect to the noise of current exper- 
iments that systematic effects are of particular concern. Consistent detection by experi- 
ments using very different techniques is important. This has been achieved only recently 
for CMB temperature anisotropy measurements, where the data obtained by DASI, CBI 

and WMAP at frequencies ^ 100 GHz are perfectly consistent with the bolometric maps 
obtained by BOOMERanG , MAXIMA, ACBAR and Archeops at 150 GHz (de Bernardis 
et al. 2003; Abroe et al. 2003; Kuo et al. 2002; Hamilton et al. 2003). 

All detections of CMB polarization to date have been made using coherent detectors 
at frequencies ^ 100 GHz. In this paper we describe a completely orthogonal experiment 
that has, for the first time, detected the CMB polarization at frequencies > 100 GHz. The 
experiment is a modification of the BOOMERanG experiment that produced the first re- 
solved images of the CMB (de Bernardis et al. 2000) and allowed the first detailed extraction 
of cosmological parameters from the CMB (Lange et al. 2001). The modified experiment, 
flown in January 2003 and hereafter referred to as BOS , is sensitive to polarization in three 
bands centered at 145, 245 and 345 GHz. We present here the measurement method, the 
instrument, and the maps of the Stokes parameters I, Q, U of the CMB detected by B03 in 
the 2003 campaign. 

Maps of CMB anisotropy and polarization are an important step in compressing the 
cosmological information into power spectra, but they are also important on their own. 
Maps are essential for understanding systematic effects in the measurement and the level of 
foreground contamination, and can be used to test the Gaussianity of the CMB fluctuations 
(see e.g. Polenta et al. (2002); De Troia et al. (2003); Komatsu et al. (2003); Aliaga et al. 
(2003); Savage et al. (2004)). 

Estimates of the power spectra (TT) , (TE) and (EE) from BOS are described in three 
companion papers (Jones et al. 2005; Piacentini et al. 2005; Montroy et al. 2005), and the 
resulting constraints on cosmological parameters in a further paper (MacTavish et al. 2005). 

2. The BOOMERanG -03 Instrument 
2.1. Generality 

This instrument derives directly from the BOOMERanG payload flown in 1997 (Pi- 
acentini et al. 2002) and in 1998 (Grill et al. 2003). That instrument provided the first 
high signal-to-noise maps of the CMB anisotropy with sub-horizon resolution (de Bernardis 



-6- 



et al. 2000; Netterfield et al. 2002; Ruhl et al. 2003), and identified three peaks in the angular 
power spectrum of the CMB (de Bernardis ct al. 2002; Ruhl et al. 2003). After the 1998/1999 
flight, the instrument was recovered and modified to make it sensitive to polarization and 
to improve the attitude reconstruction hardware. In this section we describe the different 
subsystems, with focus on the new ones. 

BOS is a scanning polarimctcr, composed of an off-axis, 1.3m diameter mm- wave tele- 
scope, a cryogenic multi-band bolometric receiver, and an attitude control system. The latter 
is able to control the azimuth and elevation of the telescope while the payload is floating in 
the stratosphere, at an altitude of ^ 30 km, under a long duration stratospheric balloon. 

We use the sky scan to modulate the signal. We map the anisotropy of the linear 
polarization by means of two separate bolometers, Bi and B2, that observe the sky through 
the same feed structure but are sensitive to orthogonal polarization directions. This device 
is called a Polarization Sensitive Bolometer (PSB). Each bolometer signal is processed and 
amplified separately. 

In principle, we can then difference the two signals to obtain the Stokes parameter Q of 
linear polarization. The U parameter is measured by means of an identical PSB, containing 
bolometers and i?4, rotated by 7r/4 in the focal plane with respect to the first one. In 
this minimal set of four bolometers, the principal axis of each sensor is rotated with respect 
to the focal plane by an angle ctfe. For the first PSB ai — and a2 — n/2; for the second 
one — 7r/4 and — Stt/A. 

In practice, our sky scan strategy uses repeated scans over the same sky pixel p. At 
different times ti during the survey, the focal plane rotates with respect to the sky by an 
angle 7^. Information on Q and U in each sky pixel thus comes from all the bolometers 
present in the focal plane, according to the relation 

K,k^\^k[Ip + QpCos[2{ak + ji)] + UpSm[2{ak + -fi)]\+ni^k ■ (1) 

Here V^^^ is the signal measured by bolometer k at time t,; is the responsivity of bolometer 
k; Ip, Qp, Up are the Stokes parameters of pixel p in the chosen celestial coordinates, and rii^k 
is the noise contribution to the i-th measurement on that pixel. This system of equations 
can be inverted and the Stokes parameters estimated if a sufficient number of measurements 
over a range of angles ctfe -|- 7^ is taken. 

Other ways to modulate the polarization involve the use of a modulating analyzer to 
extract the polarized component by synchronous demodulation. Rotating wire grids, half 
wave plates, K-mirrors, Faraday rotators, Fresnel rombs, have been used or proposed as 
polarization analyzers (e.g. Keating et al. 2002; Battistelli et al. 2002; Hanany et al. 2003; 
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Gervasi et al. 2003; Gundersen et al. 2003; Catalano et al. 2004). While all of these techniques 
can in principle provide a valuable means of reducing requirements on the stability of detector 
gains and offsets, they come at a cost in both complexity and bandwidth. Correlation 
polarimeters have so far been implemented only with coherent detectors (e.g. Carretti et al. 
2001; Padin et al. 2002). 

The polarization measurement strategy defined by eq.(l) is prone to leakage of the 
unpolarized component / into the polarized ones Q and [/, if the responsivities Sk are not 
known exactly. Similarly, errors in the principal axes angles mix Q and U into each other. 
The polarimetric calibration consists of measuring all Sk and Uk- The precision required to 
obtain the common-mode rejection needed in our case can be estimated as follows. 

For simplicity we choose the reference frame to have 7 = 0, and we consider the pair of 
detectors with ai = Q,a2 = 'n /2. In this case we can't recover all the parameters but only / 
and Q. Equation (1) becomes 
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where we have expressed the calibration constants in terms of an absolute calibration S 
and a relative cahbration TZ — S2/S1. The solution is 
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Uncertainties on the calibration constant us and on the relative calibration (Th are propagated 
in the error on Q by 
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For £'-mode polarization we expect to have a factor IjQ oi the order of 20. If we want to 
achieve a ~ 10% accuracy in the determination of Q and U, we need to have a-n/TZ ^ 1% for 
each detector. We show below that relative calibration constants can be measured within 
an error of ^ 2% by comparing the CMB temperature anisotropy measured in different 
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detectors (see §7.2.2). The absolute calibration S has to be determined with a relative error 
as/S < 5%. 

Prom equation (1) we can also estimate the acceptable principal axis angle uncertainty. 
Again we choose the reference frame to have 7 = 0, and we consider the couple of detectors 
with ai — 0,a2 — 1^/2, with errors and Uaj. To first order the measured Stokes parameter 
Qm will he Qm = Q + 2U i-e. {Qm -Q)/U ~ '^^/o%'+~o%■ From this equation 
we see that CTq. < 2° produces an error in Q (and U) ^ 10%. Mixing Q and U also mixes 
E and B-mode signals, so errors in alpha may affect the level at which one can set an upper 
limit on the B-mode polarization anisotropy signal. 

Non-ideal polarized detectors have a residual sensitivity to the polarization component 
orthogonal to the principal axis (cross-polar response). As we rotate a perfect polarizing 
grid in front of a non ideal polarization sensitive detector, the detector response is given by 
a modified Malus law: 

V = SIq{1 - {sin^(adet - agrid)), (8) 

where S is the responsivity, Jo is the incident power and ^ is the polarization efficiency, a^et 
is the angle of the main axis of the detector with respect to the x axis of the reference frame, 
and (Xgrid is the position of the principal axis of the polarizing grid. 

This can be rewritten in terms of the cross-polar response coefficient e = V±/V]| defining 
the response of the detector to radiation polarized orthogonally to its main axis, as a fraction 
of the response to radiation polarized parallel to its main axis. We have e = 1 — ^. Using an 
operative definition of responsivity <S' — S{1 + e)/2 we rewrite 



V = S'lo 



1 — e 

1 + C0S{2{adet - Olgrid)) 



(9) 



Taking into account cross-polar response, and dropping the ' for simplicity, equation (1) 
becomes 



I + Q- cos(2q;) + U- sin(2, 
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(10) 



This system can be inverted to measure Q, U, and / if the cross-polar response coeffi- 
cients are known. In the same simplified case considered above the system becomes 
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Assuming that ei and 62 are known, this is solved by 



(l + ei + ez) 



(13) 



(14) 



where the terms of the second order in e have been neglected. Proceeding as above, we find 



We conclude that in our case the uncertainty on the cross-polar response parameters can be 
^ 5%. All the conclusions obtained in this section have been verified by means of numerical 
simulations with realistic parameters (Masi et al. 2002). 

Scanning the sky with AC-coupled detectors is an effective way to measure the intensity 
and polarization of the sky at the angular scales of interest for the CMB. This strategy 
has been used by BOOMERanG -98 and Archeops, is in use in QuAD, and will be used by 
Planck-HFl and LFI. 

In the case of BOS , we scan the sky rotating the full payload in azimuth, so that 
atmospheric emission is almost constant along the scan. Working with an AC coupled 
amplifier, our system is insensitive to constant signals. 

We alternate forward and reverse scans to map the low foreground region already 
observed by BOOMERanG , centered in the constellations Caelum and Horologium, at 
RA ~ 5.5h, dec ~ —45" {b ^ — 30°). During Antarctic LDB flights, the average lati- 
tude of the payload is ~ —78°. In one day, due to sky rotation, this procedure produces a 
highly cross-linked scan pattern (see Crill et al. (2003), Fig. 9), which is important for map 
making. 



Our sky coverage is optimized to reduce the errors on the CMB power spectra given the 
constraints imposed by our telescope's hardware and the presence of bright celestial sources 
(the Sun and Galaxy). With the assumption of uniform sky coverage and uncorrelated 
noise from pixel to pixel and between T,Q and U, the errors on {EE) are approximated by 
(Zaldarriaga & Seljak 1997) 
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2.2. Observation strategy optimization 
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and the errors on {TE) are given by 
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where Nt and Ne are the effective noise equivalent temperatures in T and E of the combined 
set of detectors; r is the integration time spent uniformly covering the /sky fraction of the 
sky, and i?| is the beam response in multipole space. For the case of B03 , N"^ = 2xA^|. 
since there are eight linear detectors arranged in four orthogonal pairs, each with a different 
orientation on the sky (see Figure 5). If we use an estimate of our detector sensitivity and 
flight time (both based on the previous flights of BOOMERanG ), and a set of power spectra 
based on the best models to date, then we can calculate the uniform sky coverage needed 
to maximize the signal to noise ratio near the first peaks of (EE) and (TE) . The signal to 
noise ratios for those band-powers are a fairly fiat function of the sky coverage. The signal 
to noise ratio for (EE) peaks near Airfsky ~ 70 square degrees. For (TE) , the peak is near 
^T^fsky ~ 1600 square degrees (see Figures 1 and 2). 

The BOOMERanG telescope has a useable elevation range of 35 to 55 degrees. The 
telescope is designed to scan no more than 60° from the anti-sun direction. Exceeding this 
range could cause heating of the telescope baffies by sunlight. There are also constraints on 
the acceptable scan periods of the telescope. From the 1998 fiight of BOOMERanG , we know 
that certain scan periods excite pendulations in the balloon-gondola system. The scan speed 
is restricted by the thermal time constants of the detectors, the mechanics of the telescope 
control systems, and the stability of our readout electronics. 

Given these constraints, we created a scan strategy that came as close as possible to 
producing uniform coverage over both a "deep" region (for sensitivity to {EE) ) and a 
larger "shallow" region (for sensitivity to {TE) ) . Because each change in pointing elevation 
perturbs the telescope, we decided to adjust the elevation no more than once per hour. With 
this restriction, we found that the smallest reasonable size that we could achieve for the 
"deep" region was ~ 100 square degrees. The size of our "shallow" region {k, 800 square 
degrees) was bounded on one side by the galaxy and on the other side by the distance from 
the anti-sun direction. 

When determining the details of the scan strategy, we simulated the scanning of the 
telescope based on the same "schedule file" used to actually control the telescope during 
the fiight. These simulations produced a coverage map for a given schedule file. Since this 
coverage was non-uniform, we approximated the spectral errors by the sum of the error 



- 11 - 



3.0 



2.5 - 



^=350 



2.0 



1.5 



1.0 



0.5 



NIOOO 



0.0 







20 



40 



sky 



60 

(deg^) 



80 



100 



Fig. 1. — The signal to noise ratio vs. sky coverage for the (,EE) computed using the best 
estimate of Ne- The three curves refer to different A£ = 50 bandpowers, and are labelled 
using the center multipole. These multipoles correspond to peaks in the "concordance" {EE) 
power spectrum. 
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Fig. 2. — The signal to noise ratio vs. sky coverage for (TE) computed using the best 
estimate of Nt and A'^^. The three curves refer to different A£ = 50 bandpowers, and are 
labelled using the center multipole. The selected multipoles correspond to extremals in the 
"concordance" (TE) power spectrum. 
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contributions from each pixel 

(18) 

^x,e - ^ (2i+l)U P'^ + V ^'^ ^ r ) r^'^ + r ) 

where fp is the sky fraction of a typical pixel and Up^x and Up^E are the noise in a given 
pixel computed from the NET of the detectors and the integration time on the pixel. We 
decided to spend the first four days of the flight scanning over the "shallow" region, and to 
spend the remainder of the flight on the "deep" region. After confirming that the scanning 
schedule produced the desired sky coverage (and signal to noise), we adjusted the scan speed 
so that it did not coincide with any pendulation modes discovered in B98 . Because of these 
adjustments, B03 uses a range of many different azimuth scan speeds: between 0.2 and 0.6 
deg/s for the deep survey, and between 0.6 and 1 deg/s in the shallow survey. 

Twice a day, when the scans over the deep and shallow regions would be at nearly 
constant declination and therefore give minimal cross-linking, the telescope was directed to 
scan a region spanning the galactic plane, including RCW38 and other galactic sources useful 
for calibration (see section 7.4). 

2.3. Detectors 

Cryogenic bolometers are the most sensitive detectors for continuum mm- wave radiation 
(see e.g. Richards (1994)). 

The focal plane of B03 consists of eight optically active bolometric receivers, two dark 
bolometers, and a fixed resistor, all of which operate from a 270 mK base temperature. 
The optically inactive channels are read out using identical electronics as a check against 
microphonics, RFI and baseplate temperature fiuctuations. The focal plane is split equally 
between two types of receivers: four (dual-polarized) Polarization Sensitive Bolometer (PSB) 
pixels operating at 145 GHz, and four (single-polarization) two-color photometers using 
spider-web bolometers operating at 245 and 345 GHz. 

Polarization Sensitive Bolometers consist of a pair of co-located silicon nitride micromesh 
absorbers which couple aniostropically to linearly polarized radiation through a corrugated 
waveguide structure (Jones et al. 2003; Jones 2005). The system allows simultaneous back- 
ground limited measurements of the Stokes I and Q parameters over ^ 30% bandwidths. 
The absorbers, separated by ~ 60//m, are electrically and thermally isolated from one an- 
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other. The devices used in B03 are shown in Fig. 3. The common thermal and radiative 
environment resulting from the physical proximity of the two detectors provides gain stabil- 
ity and rejection of differential temperature fluctuations which are limited only by differences 
in the properties of the NTD Gc thermistors and the electrical leads, which determine the 
thermal conductance to the bath. Both linear polarizations propagate through a shared 
waveguide structure and set of optical filters and couple to the telescope through a cryo- 
genic corrugated feed, ensuring an identical electromagnetic bandpass and highly symmetric 
beams. PSBs are fabricated using the same proven photolithographic techniques used to 
make spider-web bolometers (Yun et al. 2003), and enjoy the same benefits of reduced heat 
capacity, low cross section to cosmic rays, and reduced susceptibility to microphonic response 
relative to monolithic bolometers. 

The two-color photometer is an evolutionary development of the photometers originally 
designed for MAX (Fischer et al. 1992), and used subsequently by the SuZIE (Holzapfel 
et al. 1997), the FIRP instrument on the IRTS (Lange et al. 1994), and the BOOMERanG 98 
(Piacentini et al. 2002; Grill et al. 2003) GMB experiments. The B03 photometer has been 
optimized for only two frequencies, and has been made polarization sensitive by the fitting 
of a polarizing grid to the feed aperture. The detectors are all similar, if not identical, to 
the detectors flown on BOOMERanG 98. The B03 feed design, consisting of a multi-mode 
back-to-back profiled corrugated horn, is significantly advanced relative to earlier versions of 
the photometer. This system achieved high efficiencies and symmetric beam patterns over 
the Ml 200-420 GHz bandwidth. 

The radiation is coupled from the photometer feed through a 420 GHz metal mesh low- 
pass filter into the 12.7 mm diameter photometer body. A dichroic filter, oriented at 45° with 
respect to the optical axis, directs radiation at frequencies above 295 GHz to the 345 GHz 
detector module while passing the lower frequencies to the 245 GHz detector module. The 
detector modules are thermally isolated from the photometer body, which is held near 2K, 
by a ~ 5-mm gap. The photometer sub-Kelvin feeds arc smooth walled with an exit aperture 
matched to the geometric area of the absorber. Corrugated feeds are not necessary, as the 
polarization discrimination and beam forming is determined by the 2-Kelvin feed antenna. 

The two configurations used in the focal plane are presented in Fig. 4. The performance 
of the receivers as integrated in the BOOMERanG focal plane is reported in Section 3. 



I Electrical Leads 





Fig. 3. — Photograph of the radiation absorber of a Polarization Sensitive Bolometer used 
in BOS . The Si^N'^ micromesh absorbing grid has a diameter of 2.6mm, the absorber leg 
spacing is 108 fim, and each leg is 3 fim wide. Only the vertical wires are metalized. The 
horizontal wires are used for structural support and do not absorb radiation. 
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Fig. 4. — Schematic of the 2-color photometers (bottom) and of the PSB polarimeters (top) 
used in BOS . 



-17- 



2.4. Optical filtering 

Optical filtering is of critical importance to a bolometric receiver; in addition to defining 
the optical pass-band, care must be taken to ensure that the detector is shielded from out-of- 
band radiation originating from within the cryostat. These filters also play a significant role 
in the determination of the end-to-end optical efficiency of the system. Finally, the optical 
filtering must reduce the radiative loading on the various stages of the cryogenic system to 
acceptable levels. 

The window of the BOS cryostat consists of a 50 — /xm polypropylene film stretched on an 
elliptical aluminum frame. The window has excellent transmission at all three wavelengths 
but is exceedingly fragile, requiring replacement after each cycle the cryostat. 

Most of the filters used in the BOS optical system consist of layers of patterned meshes 
deposited on polypropylene substrates (Lee et al. 1996), with the gaps between layers filled 
by polypropylene as well. The dichroic beam-splitter used in the photometer body is the 
sole exception, being an air-gap the inductive layers are deposited on a thin Mylar substrate 
and stretched on an aluminum frame. The polarizing grids used on the photometer and in 
laboratory testing, and the neutral-density filter, are made in a similar fashion. Instead of 
inductive or capacitive grids, a linear pattern is used for the polarizer, while a uniform lossy 
coating is used for the NDF. 

The layers of the filter are hot-pressed to form a single self-supporting filter. Some of 
the thicker hot-pressed filters are antirefiection coated with a tuned layer of PTFE (Tefion). 
PTFE has high infrared emissivity, which initially resulted in excessive heating of the BOS fil- 
ters and a large thermal load on the LN2 and "^He stages. In addition, the heating of the 
filters led to a significant increase in background loading of the detectors resulting from 
in-band thermal emission. To ameliorate this problem, large-format, composite IR blockers 
were fitted in front of the 77K and 2K filters. These filters have high in-band transmission 
and reflect radiation at wavelengths shortward of ~ lOO/zm. 

2.5. The focal plane 

To measure polarization, we combine information from spatially separated pixels, as 
shown in eq.(l). The focal plane layout (Fig. 5) is designed to minimize spatial separation 
between pixels, given the constraints of the existing BOOMERanG optics and the size of 
the feed-horns. This allows for maximal overlap of maps made by spatially separate pixels. 
The wide focal plane of the BOS telescope is thus populated by 8 pixels with independent 
corrugated feed horn systems. Each pixel contains two detectors. The four pairs of 145 
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360 GHz LPF 


410 GHz LPF 




255 GHz LPF 








540 GHz LPF 








168 GHz LPF 







^" The neutral density filter can be mechanically rotated in and out of the beam 

while the system is cooled down. 
^ A black-polyethilene (BP) Yoshinaga filter is an absorptive filter consisting of a thaUium 

salt deposited on a black polyethylene substrate, the thickness of which is 

tuned to minimize reflections (Yamada et al. 1962). 



Table 1: The optical filtering scheme employed by B03 . In order to take advantage of the low 

backgrounds available at float altitudes, much care must be taken to reduce the background orig- 
inating from within the cryostat. While the metal mesh filters, which consist of bonded layers of 
polyethylene, exhibit in-band emissivities at the percent level, the PTFE antireflection coating is 
several times more emissive. It is crucial that these filters remain well heat-sunk and protected 
from infrared emission from the warmer stages. 
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GHz PSBs in the lower row of detectors provide the best sensitivity for CMB temperature 
and polarization anisotropy. We have introduced some level of redundancy by using four 
independent PSB pairs, covering with their principal axes the range a = [0,7r] in n/S steps. 
The four pixels in the upper row are 2-color photometers operating at 245 and 345 GHz. 
These are included to provide a lever arm for discriminating CMB from dusty foregrounds. 
Table 2 summarizes the properties of the BOS receiver. 

2.6. Telescope and optics polarization properties 

The BOOMERANG telescope is an off-axis system, minimizing the radiative loading on 
the detectors. Its configuration is close to the Dragonc condition (Dragone et al. 1974, 
1982), which nulls the cross-polar response in the center of the focal plane. However, the 
need to accommodate a large number of detectors in the focal plane drove our optimization 
towards a wide corrected focal plane rather than nulling the cross-polar response only in 
the center. We optimized the optics for diffraction limited performance at 1 mm over a 
2° X 5° field of view. Radiation from the sky is refiected by the parabolic primary mirror 
(1.3 m diameter, f=1280 mm, 45° off-axis) and enters the cryostat through a thin (50 /xm) 
polypropylene window near the prime focus. Inside the cryostat, at 2K, the fast off-axis 
secondary (elliptical) and tertiary (parabolic) mirrors re-image the prime focus onto the 
detector focal plane. They are also configured to form an image of the primary mirror at the 
10 cm diameter tertiary mirror, which is the Lyot-stop of the system. In the center of the 
tertiary mirror a ~ 1 cm diameter hole hosts a thermal calibration source (callamp) which 
is flashed at flxed intervals during the flight (see Grill et al. (2003) for details). The size of 
the tertiary mirror therefore limits the illumination pattern on the primary mirror, which is 
under-flUed by 50% in area (85 cm in diameter) to improve the rejection of side-lobes. This is 
further improved by cold absorbing baffles surrounding the cold mirrors and rejecting stray 
light. Detailed parameters of the optics are described in Piacentini et al. (2002) and Grill 
et al. (2003). 

We have studied the polarization properties of this system by means of the physical 
optics code BMAX (Jones 2005). The total power beams B[9, (f) resulting from the system 
of the feed-horns and telescope are presented in Fig. 6. 

In Fig. 7 we compare the cross-polar beam (contours) to the co-polar beam (colors) 
computed for one of the 145 GHz channels. When integrated over 47r, the cross-polar response 
is ~ a few xlO~^ of the co-polar one. We conclude that the cross-polar contribution due to 
the optics is negligible with respect to the one intrinsic to the detectors. 
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Fig. 5. — Focal Plane Schematic. 2-color photometers with band centers at 245 GHz and 
345 GHz populate the upper row. Each photometer is only sensitive to one polarization. 
The lower row has 4 pairs of PSB's. The elements in a PSB pair are sensitive to orthogonal 
polarizations. The circles representing the pixels show relative beams sizes: ~ 7' for both 
photometer channels and 9.5' for the PSB's. The arrows through the circles show the orien- 
tation of the principal axis of polarization. The photometer and PSB rows are separated by 
30' in elevation, while the pixels in a row are separated by 30' in cross-elevation. The labels 
of the two bolometers used for each pixel are also reported. 
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Freq 


Bandwidth 


^detectors 


Beam FWHM 


NETcMB 


145 GHz 


45 GHz 


8 


9.95' 


170 ^Ky/^ 


245 GHz 


80 GHz 


4 


6.22' 


320 fIK^/^ 


345 GHz 


100 GHz 


4 


6.90' 


450 fiK^ 



Table 2: Summary of the properties of the BOS receiver. The noise reported in the last column is 
for a frequency of 1 Hz and is the average noise of all the detectors at that frequency. 
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Fig. 6. — Total power beams B{6, 0) computed for the BOS optical system using BMAX. The 
lower row is populated by the 145 GHz PSBs, the upper row is populated by the 245/345 
GHz photometers. Contours are plotted for power rejection levels of 0, -3, -10, -20 dB. 
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Fig. 7. — Comparison of the cross-polar (contours) and co-polar (colors) beams for one of 
the 145 GHz channels, as computed with the physical optics code BMAX. 
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2.7. Readout and signal processing 

The bolometer readout unit is the same used in the previous flights of BOOMERanG 
(Piacentini et al. 2002; Grill et al. 2003). The BOS bolometers are AC biased with a differential 
sine wave at ~ 140 Hz. The AC voltage across the bolometer is modulated by the resistance 
variations induced by changes in the absorbed microwave power. A matched pair of low 
noise J-FETs (based on Infrared Laboratories cryogenic JFET modules) inside the cryostat 
reduces the signal impedance from ~ lOMQ down to ~ IkQ; the signal is then amphfied 
by a differential preamp (AD624), band-pass filtered to remove noise outside the signal 
bandwidth, and synchronously demodulated by a phase sensitive detector (AD630). The 
output of the AD630 is proportional to the instantaneous resistance of the bolometer. High 
frequencies (above 20 Hz, i.e. above the cutoff frequency of the bolometer) are removed 
by means of a 4-pole order low-pass filter. Signal (and noise) components below 5.6 mHz 
are attenuated using a single pole high pass filter. A further 100 x amplification is applied 
to fill the dynamic range of the ADC. The total amplification of the readout chain (from 
the bolometer to the ADC input) is 50000. This ensures that bolometer plus readout noise 
is more than the quantization noise of the ADC. The output signal is analog to digital 
converted with 16 bit resolution at the sampling frequency of 60.0 Hz. The warm readout 
circuit has a gain stability of < 10 ppm/°C. The readout noise is fiat down to a few mHz, 
and is negligible with respect to the bolometer noise with NEP ~ 1.5 x 10~^^ W/\/Hz. 



2.8. Cryogenic System 

The cryogenic system is the same flown in 1997 and 1998. The main cryostat (Masi et al. 
1999) has a 65 liter tank filled with pressurized (1 atmosphere) liquid nitrogen, and a 60 liter 
liquid helium tank, which is pumped during flight so the He is superfluid. It maintains at 
1.6K a large (60 liters) experimental space containing the optical fllters, a movable neutral 
density fllter, and the ^He refrigerator (Masi et al. 1998) cooling the bolometers and feed 
horn systems. The hold time of the system is around 20 days. While for the 1998 flight 
we used two separate entrance windows for the left and right sides of the focal plane, for 
the B03 flight we used a single vacuum window, roughly elliptical (100 mm x 65 mm). The 
material is the same 50//m thick polypropylene used in 1997 and 1998. 
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2.9. Attitude Control System 

Azimuthal attitude control is provided by a reaction wheel with a moment of inertia 
approximately 0.2% of that of the complete gondola, and by a second torque motor in the 
pivot which couples the gondola to the flight train. The torque applied to the reaction 
wheel is proportional to the error in angular velocity, and the torque applied to the pivot is 
proportional to the rotation rate of the reaction wheel. On short time scales torque is from 
the reaction wheel motor, while on long time scales torque is from the pivot motor. During 
a 1 dps scan, the reaction wheel has a peak rotation rate of ~ 100 rpm. The motors are 
driven by a custom MOSFET bridge, controlled by a 20KHz PWM signal from the attitude 
control computer. The back-emf of the reaction wheel motor is compensated for in the motor 
control software, based on the reaction wheel rotation rate. 

Angular velocity readout is provided by KVH ECore2000 fibre optic rate gyroscopes, 
which provide an angular random walk noise of around 8'/v^ down to O.OlHz. Coarse 
absolute pointing is provided by a TANS-VECTOR differential GPS array (well calibrated, 
but with 6' drifts on ~ 20 min time scales), and a fixed sun sensor (sub-arc- minute precision, 
but difficult to calibrate). The position of the inner frame relative to the outer frame is 
determined using a 16 bit absolute encoder. Of these sensors, only the GPS array provides a 
complete measure of the gondola orientation in Az, Pitch, and Roll, though with significant 
uncertainty. 

The gondola is scanned in azimuth with a rounded saw-tooth wave form. At all times, 
feedback is to angular velocity from the azimuthal fibre optic rate gyroscope. During the 
linear part of the scan, the request velocity is constant (0.3°/s typical). When the end of 
the linear scan is reached (as determined by an absolute sensor, typically the differential 
GPS) the control changes from fixed angular velocity to fixed angular acceleration, until the 
velocity has reached the desired value in the opposite direction. The absolute sensors are 
only used to define the turnarounds. 

The elevation is changed by moving the inner frame, which contains the cryostat, optics, 
and receiver readout electronics, relative to the outer frame using a geared dc motor driving a 
worm gear. The elevation is only changed between observing modes. During an observation, 
the gondola is only scanned in azimuth. No attempt is made to remove pitch motions of the 
outer frame by controlling the inner frame. 

Two redundant, watchdog switched 80386 based computers take care of the attitude 
control logic, including adherence to the observation schedule file and in-fiight commanding. 
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2.10. Attitude Reconstruction Sensors 

Two pointed sensors, a pointed sun sensor and a pointed star camera, which have 
excellent intrinsic calibration, were added for the B03 flight. 

The BOS tracking star camera (SC) consists of a video, COHU brand (4920 series), 
monochrome, Peltier cooled, CCD camera equipped with a Maksutov 500 mm focal length, 
f/5.6 tclephoto lens. This setup yields approximately 4 arcsecond per pixel resolution and 
~ 30 arcminute field of view. Affixed to the lens is a seven ring baffie which is covered in 
aluminized mylar. The interior of the baffle is painted with black water-based theatre paint 
to prevent light scatter from entering the optics. For the flight, a 715 nm high pass filter 
was installed and the camera gain was set to minimum in an effort to lower the risk of CCD 
saturation at float, in the daytime Antarctic sky. 

The SC is attached to a yoke type equatorial mount. Motion control of the two axis 
system is provided by two Applied Motion Products high torque stepper motors. An encoder 
on each axis provides position feedback for controlling the motion of the mount. Motor 
current is controlled via pulse width modulation and a PID loop is used for the logic. Addi- 
tional feedback from the azimuth gyroscope is required to facilitate star tracking while the 
telescope is scanning. Video images are captured with a MATROX METEOR frame grabber 
at a rate of 10 Hz. The SC raw data consists of readouts from the two encoders, star pixel 
location in the camera field of view and star ID. Star azimuth and elevation relative to the 
gondola are reconstructed post- flight. 

The first solar sensor (Pointed Sun Sensor, or PSS) is a motorized two axis sensor (Romeo 
et al. 2002) based on a four quadrant photo-diode. Unbalance on the sensor activates the 
motors to keep the sun spot centered on the sensor, so that the sensor accurately tracks the 
sun. The angles of the sensor with respect to the payload frame are measured by means of 
two 16 bit absolute encoders (resolution 0.33'). The second Solar sensor (Fixed Sun Sensor, 
or FSS) has no moving parts (Romeo et al. 2002) and is composed of two orthogonal digital 
meridians. In such a device a hnear slit is orthogonal to a linear CCD. The sunlight entering 
the slit excites different pixels for different angles of incidence. The output of the sensor is 
the position of the center of mass of light, which can be related to the sensor-sun angle after 
appropriate calibration. Due to the variation of the luminosity level, an in-fiight calibration 
is needed. 
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3. Pre-flight calibration 

We extensively tested the instrument before the flight. This allowed tuning of some of 
the parameters for maximum performance, such as the bolometer bias frequency, preamplifier 
gains, and ACS sensor gains. It also allowed measurement of instrument parameters which 
are not expected to change from the lab to the flight, and those that cannot be measured 
in flight. These are the time-domain transfer function, the spectral response, the angular 
response, the principal axes of the polarimeters, and the cross-polar response. In the following 
subsections we describe the set of measurements we performed before launch. 



3.1. Transfer function 

In a scanning instrument, multipoles of the CMB fluctuations are encoded at frequencies 
/ in the time-ordered detector data: 

f^-i (20) 

TT 

where 9 is the sky scan speed. The frequency domain response of the instrument has to be 
known in order to deconvolve the signals obtained during the fiight. In BOS we use azimuth 
scan speeds between 0.2 and 0.5 deg/s: multipoles in the range 100 < £ < 1500 produce 
bolometer signals at frequencies in the range 0.07 Hz ^ f ^ 3 Hz. The bolometers used in 
B03 have thermal time constants in the ~ 50 ms range. Their response is closely described 
by a first order low-pass filter. The readout electronics has a transfer function that matches 
this response and cuts low frequencies in order to remove the DC level and low frequency 
drifts. The transfer functions of the detectors for the full system (bolometers -|- readout) are 
measured with a 77K blackbody load, with a cold (2K) Neutral Density Filter in the optical 
path (transmission ~ 1.5 %). This closely simulates the in- fiight loading on the bolometers, 
as confirmed by the measurement of the in-fiight DC level across the detectors. The source 
is a modulated signal obtained chopping a small 300K source against the 77K blackbody 
background (filling about 5% of the beam solid angle). The results are reported in Fig. 
8. The high frequency response is limited by the time constant of the bolometers. This is 
between 50 and 100 ms for the PSBs, and between 5 and 20 ms for the 245 and 345 GHz 
bolometers (see table 3). We estimate a ^ 10% uncertainty in applying these data to fiight 
conditions. 
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Fig. 8. — Measured transfer function of all the B03 channels. The CMB anisotropy and 
polarization signals of interest produce detector signals in the range 0.1 Hz / ^ 5 Hz. 
The low frequency response is due to the signal electronics chain. The high frequency 
response is affected by the individual bolometers' time constants; hence the variations at 
high frequencies from channel to channel. 
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3.2. Spectral Calibration 

The optical power absorbed by the bolometer is a function of the incident optical power, 
the spectral response of the bolometer and the optical efficiency of entire system. 

For an unpolarized beam filling source and a detector sensing a single polarization, the 
absorbed optical power is 

W = ^An J B{u)e{u)du, (21) 

where A is the aperture area, Q is the beam solid angle, e(i/) is the spectral response of the 
system (normalized to 1 in the maximum), 77 is the spectral normalization (i.e. r]e{i') is the 
actual spectral response), and S(i/) is the spectrum of the input source. 

The optical efficiency is then defined as the average of the function r)e{i') over the 
bandwidth of the detector. 

The spectral response of our detectors is measured using a polarized Martin-Puplett 
Fourier Transform Spectrometer (FTS). Figure 9 shows the details of the spectral response 
e{u) of the different bolometers, as measured before the flight. 

To normalize the spectral response, we use load curve data taken with the bolometers 
viewing blackbody optical loads of 77 K (liquid nitrogen) and 90 K (liquid oxygen) (Montroy 
2003). For load curves taken with two different optical loadings, we can assume that the 
incident power (optical plus electrical) on the bolometer is the same when the bolometer 
resistance measured from the two curves is equal 

Pelecihiasl) + W'"'"'^ = Pelec{Ibias2) + W''"^^ (22) 

By differencing the electrical power, we can calculated the difference in optical power 

AP = Pelecihiasl) - Pelec{Ibias2) = W^^"""' - (23) 

We calculate the spectral normafization rj from this power difference. Table 3 shows the 
spectral normalization and flat band optical efficiency calculated from load curve differences. 

The spectral calibration of the detectors allow us to compute color corrections for all our 
photometric channels. In fact our system is calibrated in flight using the CMB anisotropy 
as a reference source (see §7.2.2). For this reason the responsivity must be corrected when 
observing sources with a different spectrum, like, e.g. interstellar dust emission. However, 
the spectral matching of PSBs is good enough that even not applying any color correction, 
the spurious polarization degree is ^ 5 x 10""^ for any reasonable ISD spectrum. 
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Fig. 9. — Measured spectral response e(z/) of the 16 bolometers used in BOS . The spectra 
have been measured before the flight by means of a thermal blackbody source and a Martin- 
Puplett interferometer. All the transmission measurements are normalized for unit peak 
transmission. The efficiencies of the different channels are reported in table 3. The dashed 
line is the normalized brightness of CMB fluctuations, normalized to the peak (217 GHz) 
brightness: our three bands bracket the maximum of such a spectrum. 
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Channel 


{u) (GHz) 


ViiM') (GHz) 


MJy/sr 

KoMB 


Spec. Norm. 


Opt. Eff. 


T (ms) 


145W1 


147.3 


13.5 


388 


0.35 


0.26 


89 


145W2 


146.7 


12.4 


387 


0.35 


0.26 


50 


145X1 


146.2 


12.6 


386 


0.37 


0.28 


58 


145X2 


146.6 


12.4 


387 


0.33 


0.25 


43 


145Y1 


147.0 


12.8 


388 


0.38 


0.29 


97 


145Y2 


146.8 


12.7 


387 


0.38 


0.31 


82 


145Z1 


147.0 


12.6 


388 


0.26 


0.19 


66 


145Z2 


147.1 


13.1 


388 


0.21 


0.14 


81 


245W 


248.0 


26.0 


461 


0.38 


0.33 


12 


245X 


250.7 


24.7 


460 


0.36 


0.28 


20 


245Y 


244.0 


24.6 


465 


0.33 


0.27 


15 


245Z 


247.6 


26.5 


461 


0.32 


0.27 


16 


345W 


340.3 


38.6 


321 


0.81 


0.70 


7.9 


345X 


338.9 


38.9 


323 


0.89 


0.77 


8.0 


345Y 


344.1 


37.6 


314 


0.82 


0.72 


4.7 


345Z 


337.8 


38.7 


326 


0.87 


0.68 


12 



Table 3: Spectral and time response calibrations: For all channels (named in column 1) we report 
in the second column the average frequency, and in the third column the optical bandwidth. These 
are estimated in terms of the first and second moments of the transmission spectra of Fig. 9. The 
conversion factor between Specific Brightness and CMB temperature fiuctuations, as computed from 
the transmission spectra of Fig. 9, is reported in the fourth column. Spectral normalizations and 
fiat band optical efficiencies are measured using NDF-up load curve power differences. The spectral 
normalization rj (fifth column) is calculated using the spectral response from the FTS measurements. 
The optical efficiency {r]e{u))i,and (sixth column) is calculated assuming a flat spectral response. 
The detectors are assumed to be single-moded. The high optical efficiency of the 345 GHz channels 
is likely due to propagation of multiple modes to the detector. The time constant r is reported 
in the seventh column, as measured in the laboratory with loading conditions similar to the flight 
ones. 
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3.3. Polarimetric Calibration 

To characterize the polarization properties of our system, we measured the polarization 
efficiency, ^ (or the cross-polar response e = 1 — and the polarization angle of each 
detector, a^et (see eqn. 8 and 9). 

We made these measurements with the receiver alone and with the whole system in- 
cluding the telescope. 

Figure 10 shows the configuration of the receiver polarization calibration. A rotating 
polarizing grid is placed directly beneath the cryostat. It sits above a cold load (liquid 
nitrogen) which is modulated by a chopper wheel. The chopper wheel is rotated at 2 Hz and 
the grid has a rotation period of approximately 10 minutes. When the transmission axis of 
the grid is aligned with the polarization axis of the detector, the detector sees whatever is 
behind the grid (i.e. the chopper wheel or the cold load). When the grids are 90° out of 
alignment the detector sees radiation reflected off the grid. The size of the source aperture 
is such that the modulated signal fills ~ 5% of the receiver sohd angle. The best fit values 
of the cross-polar response parameter e when the source is on the optical axis are reported 
in table 4. 

For the measurement of the full system we used a fully polarized modulated plane 
wave filling the telescope aperture. This was produced using a thermal source in the focus 
of the spare BOOMERanG primary mirror. This inverted telescope system was placed in 
front of the BOS telescope (see Figure 11). Analytical calculations (Djik et al. 1974) and 
optical simulations show that the cross-polar response introduced by the spare primary at our 
frequency is ^ 1%. The source at the focus is made with a high efficiency wire grid (12.5 fim 
diameter Tungsten wires, 25 /im spacing) in front of a chopper alternating a 77K blackbody 
and a 273K blackbody. The wire grid axis was tilted by 28° with respect to the optical axis of 
the system. The B03 telescope was mounted on a rotating platform with absolute encoders 
for angle measurements. Bolometer data and the chopper reference signals were read by the 
Data Acquisition System of the experiment. A software demodulator measuring amplitude 
and phase of the synchronous signal was used to extract the polarized source signal from the 
background. Data were taken for 18 positions of the wire grid between and 180 degrees. 
After correcting for the tilt of the grid with respect to the optical path, the data for all the 
bolometers and for all the incidence angles were ffi with equation (9). 

The measurement was carried out in two different configurations. In the first we used 
a wide aperture on the cold source, so the telescope beam was completely filled by a ^ 1° 
source. The source-telescope relative position was then optimized for maximum signal. In 
this configuration we measured the integrated polarization properties of the instrument. The 
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Fig. 10. — The source used for the polarization cahbration of the receiver. The system is 
tilted at 22.5° so that rays reflecting off the grid go to 300 K eccosorb absorbers (shown as 
black, thick lines in the drawing). The aperture of the system is ~ 5 cm. The chopper wheel 
chops between 77 K and 300 K, while a belt drive rotates the polarizing grid. The chopper 
wheel rotates at 2 Hz, while the grid has a rotation period of about 10 minutes. 
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Fig. 11. — The apparatus for measuring polarization angle and cross polar response of the 
entire telescope. The cold load is kept at 77K by a hquid nitrogen bath. The source telescope 
is a 1.3 m off axis paraboloid (the spare of the BOS mirror), and is kept at a distance of ~ 6 
m from the BOS telescope. The source is scanned in Azimuth by rotating the BOS payload 
around its vertical axis. With a small diaphragm in the source, scans are repeated at different 
elevations to obtain a full map of the beam and of its polarization properties. With a large 
diaphragm, producing a beam-filling beam, we change azimuth and elevation of BOS to find 
the maximum signal. 
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measured data for this configuration are shown in Figure 12 for the PSB channels, and in 
Figure 13 for the two-color photometer channels. The best fit values of the parameters e 
and adet arc reported in table 4. The measured cross-polar response for the PSB channels 
is of the order of 0.1. It is lower for the two color photometer channels, because of the high 
efficiency of the wire grids. 

In the second configuration we used a small aperture (~ 3' FWHM), in order to in- 
vestigate the cross-polar response properties as a function of the off-axis angle. In general 
the on-axis cross-polar response is lower than the integral values of table 4, while off-axis is 
higher. This is in agreement with the physical optics model of our system (shown in Fig. 7) 
and with the receiver polarization calibration described above. 

3.4. Beam Measurement 

We measured the beam profile B{6, (p) of our telescope with a thermal source tethered 
on a balloon in the far field of the telescope. The measurement was made in the relatively 
transparent atmosphere of Antarctica. A small tethered balloon lifted a microwave absorber 
(made out of cccosorb) into the telescope beam, at a distance of about 1.5 km from the 
payload. We used two sources to measure the beam, a sphere (with 45 cm diameter) to 
measure the main beam, and a larger cylinder (with 76 cm diameter and 91 cm height) to 
map the sidelobes. The telescope was steered by a combination of a slow azimuth scan with 
a very slow elevation drift to map the beam profile of all the BOS detectors. The position of 
the source relative to the telescope was continuously measured with the tracking star camera. 

As an example, we plot in Fig. 14 one of the 145 GHz beam profiles B{9), and compare it 
to the beam computed from the BMAX physical optics model of our system. The agreement 
is very good. This measurement, together with scans over point sources during the fiight, 
gives us confidence that the beam computed with BMAX (plotted in Fig. 6) is accurate for 
all our purposes. 

From the BMAX beam we compute the transfer function of the instrument in multi- 
pole space. This is plotted in Fig. 15, where it is also compared to a gaussian beam transfer 
function. This -B| has been convolved with the smearing resulting from pointing jitter (see 
§5.1), and used to deconvolve all our power spectrum measurements, as reported in (Jones 
et al. 2005; Piacentini et al. 2005; Montroy et al. 2005). 
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0.062 ± 0.005 
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0.11 ± 0.02 
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1 70 /I _|_ 9 n 


0.055 ±0.013 


0.09 ±0.02 


145X2 


90.0 


88.7 ±2.0 


0.070 ±0.008 


0.13 ±0.02 


145Y1 


157.5 


158.6 ±2.0 


0.051 ±0.012 


0.11 ±0.02 


145Y2 


67.5 


66.2 ±2.0 


0.060 ± 0.006 


0.11 ±0.02 


145Z1 


112.5 


109.8 ±2.0 


0.182 ±0.004 


0.22 ±0.02 


145Z2 


22.5 


21.2 ± 2.0 


0.088 ± 0.008 


0.12 ± 0.02 


245W 


135.0 


139.2 ±2.0 


0.007 ±0.005 


0.02 ± 0.02 


245X 


45.0 


42.9 ±2.0 


0.007 ±0.005 


0.02 ±0.02 


245Y 


180.0 


178.7 ±2.0 


0.000 ±0.015 


0.03 ±0.02 


245Z 


90.0 


85.2 ±2.0 


0.014 ±0.013 


0.01 ±0.02 


345W 


135.0 


139.9 ±2.0 


0.008 ±0.005 


0.01 ±0.02 


345X 


45.0 


42.6 ±2.0 


0.008 ±0.001 


0.01 ±0.02 


345Y 


180.0 


178.7 ±2.0 


0.004 ±0.014 


0.02 ± 0.02 


345Z 


90.0 


84.9 ± 2.0 


0.019 ±0.013 


0.02 ± 0.02 



Table 4: Cross-polarization and principal axis angle measurements for all the channels. The 
polarization angle and the beam filling values refer to measurements on the complete in- 
strument obtained with the apparatus described in Fig. 11, while the on-axis measurements 
refer to measurements on the receiver alone, obtained with the apparatus described in Fig. 
10 and also in a test cryostat testing single channels. 
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Beam- integrated polarization response of the entire instrument (PSB channels). 
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Fig. 13. — Beam-integrated polarization response of the entire instrument (two-colors pho- 
tometers.) 
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Fig. 14. — Measurement (crosses) of the azimuth beam profile of channel 145W1, obtained 
from repeated scans of a thermal source in the far field of the telescope. The measurements 
fit very well the prediction of the physical optics code BMAX (dashed line) . 
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Fig. 15. — Upper panel: Instrument response at 145 GHz in multipole space, computed from 
the BMAX physical optics model (red line). The beam is similar to a 9.8 arcmin FWHM 
gaussian (black line). The ratio between the window function for the actual beam and that 
for a gaussian beam is plotted in the lower panel. 
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4. Observations 

The instrument was launched by NASA-NSBF on Jan. 6, 2003, from Wilhams Field, 
near the McMurdo Station, in Antarctica (167° 5.760' E ; 77° 51.760' S). The flight lasted 
until Jan. 21, with a total of 311 hours. 

The altitude of the payload during the flight is reported in Fig. 16. The periodic 
variation is due to the daily change of elevation of the sun, while the long term trend is due 
to a small leak in the balloon. We released ballast, in a moderate amount on day 3, and then 
in a full drop on day 5. After day 11, the altitude dropped below 28 km, telescope pointing 
become difficult, and we had to stop observations. 

We devoted a total of 119 hours to scans on the deep survey region, a total of 79 hours 
to scans on the shallow survey region, and a total of 30 hours to scans over the Galactic 
plane. The remaining hours have not been used for the data analysis, due to spurious signals 
after events like ballast drop or elevation changes, and due to testing, cryogenic operations, 
non optimal performance of the attitude control system at the lowest altitudes. 

Maps of the sky coverage for the three surveys are reported in Fig. 17 and in Fig. 18. 
The histogram of the integration time per pixel for the 3 different regions is plotted in Fig. 
19. 



4.1. Attitude Control System performance in flight 

Motion control was flawless in both the azimuth scanning and the drive of the inner 
frame elevation when the payload was at nominal altitude ( ^ 25km). The GPS and FSS 
provided the required coarse, ~10 arcminute in-flight azimuth pointing and the gyroscopes 
functioned almost continuously, aside from short periods of dropouts or computer reboots. 
Near the end of the flight, once the payload altitude dropped below ~ 23 km, high winds 
made motion control impossible and the ACS was shut down. 

Performance of the tracking sensors, the SC and the PSS, was marred by communication 
problems in the ACS flight computer. With both sensors operating at the same time the 
computer was unable to parse the large amounts of incoming and outgoing data. To resolve 
this issue the PSS was turned off early in the flight. The SC provided the flne pointing data 
for the flrst and last thirds of the flight. The cold temperature at low (pre-ballast drop) 
altitude on day 4 caused the SC to freeze up, at which point the PSS was turned back on, 
providing flne pointing for the middle third of the flight. 
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Fig. 16. — Bottom: The altitude of the payload during the balloon-flight. The sinusoidal 
variation is due to the daily change of elevation of the sun, while the long term trend is 
due to a small leak in the balloon. The arrows indicate ballast releases used to counter the 
downward drift. Top: Drift of detector responsivity during the flight as measured with the 
on-board calibration lamp: it is contained within a few % for the whole length of the flight. 
The correlation with the altitude is due to to the variation of atmospheric pressure with 
altitude, which changes the temperature of the superfluid ^ife bath in the cryostat, with a 
consequent change in bolometer temperature. 
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Fig. 17. — Map of the observation time (integral over all the 145 GHz bolometers) in the 
" Shallow" and " Galaxy" surveys of BOS . The larger box includes pixels of the shallow survey 
actually used for the power spectrum analysis; the smaller box refers to the "Deep" survey 
(see Fig. 18 ). 
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Fig. 18. — Map of the observation time (integral over all the 145 GHz bolometers) in the 
"Deep" survey of BOS . This region overlaps completely with the central part of the shallow 
survey region. The box includes pixels of the deep survey actually used for the power 
spectrum analysis; however, the three circular regions marked around strong AGNs have 
been excised. 
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Fig. 19. — Histogram of observation time per pixel (3.5' pixels) in the three surveys performed 
during the BOS flight. 
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5. Raw data processing 

The polarization signal we are looking for is very small. In the observations of the deep 
region at 145 GHz the contribution from CMB polarization is of the order of a thousandth 
of the noise for a single signal sample, and of the order of a tenth of the error for the final 
measurement of a single 3.5' pixel. For this reason we must properly exclude instrumental 
artifacts and avoid the introduction of bias. We developed two completely independent 
analysis pipelines, which will be called IT (Italy) and NA (North America) in the following 
discussion. In this section and in the following we outline the data analysis and specify where 
the two pipelines make different choices. Despite the many differences, the final spectra and 
maps are fully consistent. This is the best test of the robustness of the dataset against 
alternative data selection and analysis procedures. 

The raw data acquired by the B03 telescope is convolved with a transfer function con- 
sisting of the thermal response of the detector and the filtering of the readout electronics. 
It also contains transients caused by cosmic rays and our calibration lamp, and longer time- 
scale events caused by thermal instabilities resulting after each elevation change. At each 
step of the data cleaning, we fiag (mark as bad) some samples in the datastream and fill 
these flagged samples with a constrained realization of the noise. This noise filling uses lin- 
ear prediction to replace flagged samples based on the unflagged data on either side. After 
removing very large cosmic rays, we flag a minimal number of samples around each elevation 
change and flt out a decaying exponential to the signal recovery. 

Next, the datastreams are deconvolved with the combined thermal and electrical transfer 
functions. For the NA pipeline, we used transfer functions determined by pre-flight lab 
measurements (described in §3.1). For the IT pipehne we used the in-flight determination 
of the transfer function described in §6.1. 

The samples containing signals from our calibration lamp are flagged. Medium ampli- 
tude cosmic rays are then flagged by passing a simple template through the timestream and 
looking for spikes larger than a certain size. The data is then heavily band-pass flltered to 
make small cosmic rays more visible. A second template is used to flag these small cosmic 
rays. This process might mistake bright galactic sources for small cosmic rays. To avoid 
this, we make a list of all bright source crossings in our datastreams by using the coincidence 
of sources in both detectors of a pixel pair. These samples are explicitly unflagged. After 
building up a list of flagged samples, these samples are filled with noise in the deconvolved 
data. 

After completing this primary data cleaning, we found that each of our time-streams 
contained an extraneous component that was well correlated with the accelerations of the 
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telescope's pitch and roll gyroscopes. We fit these accelerations to the data in hour long 
chunks and subtracted this contribution from the time-streams. 



5.1. Post-flight Attitude Determination 

Pointing reconstruction consists of determining where each detector was pointing in 
the sky at each data sample. We had a redundant number of attitude sensors that can 
be combined in different ways to obtain the final pointing solution. The two pipelines are 
significantly different in this respect. 

5.2. IT Pointing reconstruction 

Pointing is described by the rotation matrix S that transforms from the Celestial refer- 
ence frame to the telescope reference frame. In addition, the constant offsets of the pointing 
directions of the horns in the telescope frame must be determined. The problem is separated 
into two parts S = A- B, where the rotation matrix A converts from an Earth local reference 
frame to the telescope frame, and B is the astronomical conversion. A is called the attitude 
matrix and is defined by three Euler angles. It can also conveniently be described by a 
quaternion. 

To derive the attitude matrix we use several attitude sensors, measuring either absolute 
angles or angular velocities. In our case the absolute sensors are a differential GPS, a Stellar 
Sensor and two Solar Sensors. The differential GPS can by itself measure the attitude matrix 
A, but suffers for low frequency drifts, of the order of 7' / \fh. The solar and the stellar sensors 
measure the sensor-Sun or sensor-Star angles in the payload frame, each providing two 
equations for the three unknown Euler angles. Their combination thus provides a solvable 
system. The angular velocities are obtained from the laser gyroscopes, which have a bias 
due to their orientation in the Earth's magnetic field. 

We combine the signals of the different sensors taking into account their noise/drift 
properties, by means of an optimized Kalman filter (Kalman 1960) and by properly flagging 
the sensors signals. The Kalman filter is based on the propagation from the status at time t 
to the status at time t -|- At as predicted by the dynamic measurements (from gyroscopes) , 
in combination with the direct measurements of the status i -|- At (update). The weighting 
of the update is given by the noise properties of the gyroscopes and of the sensors. 

The direction of hues of sight of the different detectors in the telescope frame are derived 
from observation on known sources, i.e. AGNs and Galactic HII regions. 
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A final correction to the IT pointing solution is obtained by dividing the timestream in 
hour-long chunks. For each chunk we make a sky map, and find the Azimuth and elevation 
offsets which best fit a single template map derived from the BOOMERanG data from the 
1998 flight. An IRAS template was used for the scans on the Galactic plane. The pointing 
solution is then corrected using these best fit offsets, which are of the order of ^ 1 arcmin. 

The accuracy of the IT pointing solution is tested by means of the signals detected in 
the direction of bright AGNs (see §6.2), and is of the order of (2.4 ± 0.3)' rms. 

5.3. The NA Pointing Solution Pipeline 

Spikes in the raw pointing data are removed, and where possible (over a few samples) 
the data can be hnearly interpolated. Post-flight re-cahbrations are applied to the SC, PSS 
and FSS data. For the SC/PSS azimuth and elevation this is a simple matter of rotating 
the coordinates through a small angle until correlated signal is minimized. The FSS is re- 
calibrated using azimuth data from the GPS and the SC. A look-up table of sun elevation 
versus raw FSS azimuth is constructed. Each element of the table contains a GPS/SC derived 
sun azimuth relative to the gondola, averaged over the whole flight. Raw FSS data is replaced 
by the corresponding element in the table. Another useful derived quantity is obtained from 
the GPS up, north and east velocity data. The relatively stable LDB environment allows 
one to model the gondola as a pendulum. With this model first estimates of gondola pitch 
and roll, pitchov and rollav may be calculated from the GPS velocity data. 

Clean, calibrated pointing fields are combined to determine the gondola attitude. The 
sensors used in the final analysis include the FSS azimuth, the SC elevation, pitchov, rollav 
and the integrated azimuth, pitch and roll gyro data. At frequencies below ~50 mHz the 

pointing solution is based on the best fit azimuth, pitch and roll to the sun and star positions 
as determined by the star camera, the FSS, pitchcv and rollcv- Gaps in the pointing data 
less than ~ 40 seconds long are filled with integrated gyro data and above ~ 50 mHz the 
pointing solution is strictly gyro signal. Gaps longer than ~ 40 seconds are fiagged. 

The elevation encoder signal is added to the gondola pitch, thereby translating gondola 
(outer frame) attitude into telescope (inner frame) attitude. The beam offsets for each 
detector are obtained from fits to the five brightest QSOs in the CMB field. Galactic and 
CMB source centroid offsets reveal a 0.1° shift in gondola pitch after the mid-flight ballast 
drop. To account for this approximately 6 hours of data during and after the ballast drop are 
flagged and a pitch shift is apphed to all pointing data preceding the drop. The reconstructed 
elevation and azimuth of each beam on the sky, along with the measured polarization angles 
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for each detector, and the GPS latitude, longitude and time are combined to determine the 
right ascension, declination and psi angle for each beam. 

The CMB field pointing error, based on comparisons of analytical beams with observed 
beams, is ~2.5' rms in azimuth and ~1.5' rms in elevation (see §6.2). 

6. In-flight Calibration 

6.1. Transfer function 

The transfer function of each detector, including the readout, has been measured in 
flight using a procedure similar to the one described in Grill et al. (2003). Even though 
the detectors are designed to minimize cosmic rays cross section, each detector produces 
detectable cosmic ray hits at the rate of about one every two minutes, leaving on the data- 
stream a typical signature which is the response of the system to an impulsive input. This 
signature is the transfer function of the system in real space. 

Of all the cosmic rays events in a given channel, only the subset producing a spike in 
the (0.5-4) Volts range is selected for this analysis, in order to neglect the effects of noise 
and to avoid saturated events. Each event in the database is shifted in time and normalized 
to minimize the chi square to a first order approximation of the impulse response. The shift 
is performed on a grid much finer than the 60 Hz sampling rate of the Data Acquisition 
System. The combination of all the hits provides a template of the impulsive response of 
each detector. As an example, we plot in Figure 20 the data for detector 145 Wl. The result 
is insensitive to the choice of the first order approximation. The Fourier transforms of those 
templates arc the transfer functions of the detectors. The method is sensitive in a frequency 
range between ~ 0.1 and ~ 250 Hz. The resulting transfer function is similar (but not 
exactly equal) to the one described in §3.1. The small differences are probably due to the 
internal time constant of the detector absorber, which affects in different ways the response 
to mm-wave photons (depositing their energy on the entire absorber) and the response to 
particles (depositing their energy on a small localized part of the absorber) . The IT pipehne 
uses this transfer function, complementing it with lab calibration data (§3.1) at / < O.lHz. 
The NA pipeline uses the transfer function derived from pre-fiight measurements. 
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6.2. Effective Beam Calibration 

The effective BOS beam profiles are estimated combining physical optics modelling and 
observations of bright extragalactic point sources. 

The physical beam is computed with a physical optics simulation of the telescope (see 
§2.6 and Fig. 14); its shape is similar to a Gaussian function with a FWHM of 9.8 arcmin 
for the 145 GHz channels and of 5.6' for the 245 and 345 GHz channels, with residuals of 
the order of 2% (see Fig. 15). 

Pointing errors produce a jitter on the map that we assumed to be gaussian and uniform 
in all the observed sky regions. Under this assumption the jitter is defined by a single 
parameter, ajuter- 

The observed profile of a point source is the convolution of the physical beam with the 
gaussian describing the pointing jitter. For a gaussian and isotropic physical beam we would 
have 

^2 



B{9, (T jitter) = exp 



(24) 



2('^6eam + '^jitter). 

Using the actual physical beam gives a result very similar to equation 24. We measured 
Cjitter by comparing the expected beam to the profiles of five sources derived from the 8 
PSBs maps at 145 GHz. We did not use the higher frequency channels to estimate the jitter 
because of their higher noise. The jitter, estimated as the weighted average of the 5 best fit 
values, is CTjuter = (2.5 ± 0.3) arcmin. The FWHM of the effective beam used in the final 
analysis is (11.2 ± 0.3) arcmin for the 145 GHz channels, (7.8 ± 0.3) arcmin for the 245 and 
345 GHz channels. 



6.3. Gain drift from Calibration Lamp 

The responsivity of the bolometers depends on the operating temperature and on the 
radiative loading. Both these quantities can change during the flight. 

For this reason, we include a calibration lamp in the cold optics of BOOMERanG , 
that flashes every ~ 15 min. during the flight. Details on the lamp, which is a composite 
bolometer structure suspended in the center of a 1 cm hole in the tertiary mirror, are in Grill 
et al. (2003). Its intrinsic stability depends on the fact that the temperature achieved during 
the flashes is signiflcantly higher than the 2K base temperature. In this limit, the emission 
of the lamp depends only on the heater current, which is stable to better than 0.1%. 

The amplitude of the calibration signal changes only by a few percent during the flight. 
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as shown in the top panel of Fig. 16. To demonstrate that these fluctuations are not intrinsic 
to the source, we have analyzed the DC signal across the bolometers. This is proportional 
to the resistance of the bolometers, which in turn is proportional to the responsivity. For 
all the detectors we found a very good correlation between the bolometer responsivity esti- 
mated from the calibration lamp and the DC level: a clear indication that the signal of the 
calibration lamp is constant through the flight. 

6.4. Bolometer Performance 

We compute in-fiight bolometer parameters with a bolometer model based on Holmes 
et al. (1998), Jones (1953), and Mather (1982) (sec Jones (2005) for details). The receiver 
model computes the total optical load on the detectors using the resistance and thermal 
conductivity as functions of temperature measured on the ground, and the in-fiight DC 
bolometer voltage, bias current, and noise. 

We further model the optical load on the detectors due to the atmosphere, telescope, 
and Cosmic Microwave Background, and assume that leftover optical load is due to internal 
parasitic loading. The optical load and responsivity change by several percent during the 
flight (sec §6.3) due to changes in emission from the atmosphere and thermal radiation from 
cryostat components. 

We model the contributions to the noise from photon shot noise, Johnson noise, phonon 
shot noise in the thermal link between the bolometer and the cold stage, and amplifier noise 
and compute the NEP (Noise Equivalent Power) contributed by each. The results of these 
calculations, for typical values of the in-flight optical load, are reported in table 5. 

A source of impulsive noise is cosmic rays. All the BOS detectors have been designed to 
minimize their cross section to cosmic rays. As a result, even in the polar stratosphere, the 
average rate of cosmic-rays hits in PSBs is one every two minutes. Over all the identified 
events, 83% produce transients with amphtude at the ADC 0.05 < A < IV; 10% have 
amphtude IV < A < 2V, 5% have 2V < A < W and 3% have A > W. All these events 
are flagged and removed from the analysis as explained in §5. 

6.5. Noise Estimation 

The measurement of the actual in-flight noise is obtained via a Fourier transform of 
the raw datastream. A sample power spectrum of the raw data (deconvolved from system 
frequency response, as estimated in §3.1 and §6.1) is shown in Fig. 21. Signal contributions 
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units 




145 


245 


345 


Optical power 


pW 




0.496 


0.846 


6.158 


Equivalent R-J temperature 


K 




4 


4 


16 


CMB /total 


% 




23.0 


6.2 


0.9 


Atmos. /total 


% 




4.1 


32.6 


64.7 


Refl./total 


% 




21.4 


17.2 


10.0 


Internal/total 


% 




51.5 


44.0 


24.4 


NEPtotal 


10-^^ W/^ 


/Hi 


2.49 


3.89 


8.95 


NEPphoton 


10-1^ w/^ 




0.99 


1.64 


5.31 


Bose term 


% 




14.8 


7.2 


16.4 


NEP V, 

i N phonon 


10-1^ w/^ 




1.003 


1.636 


2.891 


■^EPjohnson 


10-1^ w/^ 


/Hi 


1.307 


1.839 


3.669 


■^EPamplifier 


10-1^ w/^ 


/Hi 


1.578 


2.518 


5.494 


NEPixjIo / NEP background ■ ■ ■ 






1.67 


1.50 


0.88 


NEED 


mJy/VHz 


37 


43 


177 


NETcmb 


A^KcMBA 




151 


341 


551 


NETrj 




s 


89 


92 


51 


NET-^, 






53 


171 


275 


Voltage noise 


nV/\/ife 


18.0 


17.7 


18.5 



Table 5: The receiver performance for a typical channel from each of the BOOMERANG bands, 
computed from a receiver model using the observed average in-flight loading conditions. The 
performance measured in flight is consistent with these estimates. 
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are present above the noise, at the scan frequency and at its harmonics, and in the ~ 0.1 
Hz range. The Generahzed Least Squares (GLS) map-making method rehes on knowledge 
of the noise correlation function, which is related to the noise power spectral density by a 
Fourier transform. This quantity must be estimated from the data themselves, which are a 
combination of noise and signal. 

The standard way to characterise the underlying noise is by subtracting an estimate 
for the signal (see Ferreira & Jaffe (2000)). Since the signal can be estimated by making a 
map, the two problems are entangled, and a possible approach is iterative (see Prunet et al. 
(2001) and Dore et al. (2001); see however Natoli et al. (2002) for an alternative approach). 

In the IT pipeline, a rough estimate of the signal is obtained by naively coadding a 
band-pass filtered version of the timeline; the resulting map is used as a baseline to obtain 
an estimate of the noise power spectral density. This is in turn used to make a new, GLS, 
map and the process can be iterated as desired. We have found empirically that, for BOS , 
no significant improvement is obtained by iterating more than five times. Furthermore, tests 
on simulated data suggest that, under these conditions, the underlying noise properties are 
recovered without any substantial bias (see de Gasperis et al. (2005)). 

The following scheme is employed to optimally use the information produced by all 
detectors. The above procedure is repeated for each uncahbrated bolometer timeline (e.g. 
eight times, one for each 145 GHz bolometer). Then, a first estimate of the relative calibration 
factors between channels at a given frequency are computed as explained in §7.2.1. A multi- 
channel, high signal to noise ratio, relatively calibrated map is produced and used as an 
estimate for the underlying signal to evaluate a more precise noise power spectral density 
for each bolometer. As a further step, relative calibration factors are recomputed with these 
noise estimates, and a new multichannel map is made. We have verified that iterating further 
is useless since if this last map is used to produce further noise estimates, the latter do not 
change significantly. 

A sample noise power spectrum is shown in Fig. 21. With this method we estimate the 
noise of all the BOS channels in V j \J Hz. This is converted into NEP (see table 6) using the 
measured in-fiight responsivity (see §7.2.1 and §7.2.2). 

There is a component of noise correlated between different detectors. This is especially 
visible at very low frequencies. In Fig. 21 we include an example. This noise has a similar 
level for all detector couples (roughly independent of whether the detectors arc members of 
the same PSB pair). It may originate in scan synchronous microphonics and in fluctuations 
of the residual atmosphere. The levels of correlated noise in the interesting frequency band 
(0.1-2 Hz), for all possible couples of 145 GHz bolometers, are compared to the levels of 
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Channel 


NET (iiKcmbVs) 


145W1 


140 


145W2 


137 


145X1 


156 


145X2 


150 


145Y1 


151 


145Y2 


159 


145Z1 


182 


145Z2 


281 


245W 


281 


245X 


358 


245Y 


316 


245Z 


331 


345W 


459 


345X 


424 


345Y 


620 


345Z 


309 



Tabic 6: In flight detectors performance: We report the noise equivalent temperature at 1 
Hz, derived from the noise equivalent voltage measured in flight as in §6.4, from the responsivity 
estimated as in §7.2.1 and from the spectral calibrations of §3.2. In our survey, this frequency 
corresponds to multipoles 400 ^ i ^ 1200. The noise increases at lower and at higher frequencies, 
as shown in Fig. 21. 
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noise of each bolometer in table 7. The presence of this correlated noise is neglected in the 
map-making, but is taken into account for our estimates of the power spectra (see Jones 
et al. (2005), Piacentini et al. (2005), Montroy et al. (2005)). 



6.6. CMB dipole 

The dipole of the CMB is visible as an approximately hnear drift along our short scans. It 
produces a scan-synchronous triangle wave (see Fig. 22) that we have removed from the time- 
streams before proceeding with the analysis. Since it is scan-synchronous, as some possible 
systematics, we cannot use it for a precise calibration of the instrument: the calibration 
accuracy achievable with the dipole is ^ 15%, not sufficient for our purposes. However, 
its amplitude is consistent with the amplitude of the CMB Dipole measured by COBE and 
WMAP . 



7. Stokes Parameters Maps 

7.1. The Map Making Pipeline 

A Generahzed Least Square (GLS) method is used to jointly estimate the Stokes param- 
eter sky maps /, Q and U from time ordered data. The best maps are obtained by combining 
data from all detectors available at a given frequency. This approach is implemented in the 
ROMA map making code (de Gasperis et al. (2005)) for the IT pipehne, and in the DIQU 
code (Jones et al. 2005b) for the NA pipehne. Following eq. 10, we assume the following 
data model: 



T)» — A^ 



Ip + - — — (Qp COS 2al + Up sin 2al) 
1 ~l~ 



+ nl. (25) 



Here i, t and p label channel, time and map pixel respectively; Dl are the time ordered data 
for a given channel, related to the sky maps Ip, Qp and Up by the pointing operator Alp, 
which we assume to be as simple as possible: its elements are equal to one when a pixel is 
observed, and are zero otherwise. Eq. (25) can be recast in a more compact formalism by 
defining a generalized pointing matrix Atp which includes the trigonometric functions and a 
map triplet Sp = {Ip, Qp, Up): 

Vt = AtpSp + n. (26) 
The GLS estimator for Sp is then (see de Gasperis et al. (2005)) 

Sp = (A*N-i A) A*N-ip, (27) 
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Fig. 20. — In-flight response of the 145W1 channel to an impulsive event. The frequency 
response of the system is the Fourier Transform of this response. The points are accumulated 
from several cosmic-rays events shifted and normalized to fit the same template. 



Tabic 7. 



Wl W2 XI X2 Yl Y2 Zl Z2 



Wl 


1.0 














W2 


6.0E-02 


0.85 












XI 


3.4E-02 


3.1E-02 


1.1 










X2 


3.3E-02 


3.1E-02 


5.9E-02 


1.0 








Yl 


4.8E-02 


4.0E-02 


3.5E-02 


3.5E-02 


1.2 






Y2 


4.3E-02 


3.6E-02 


3.5E-02 


3.5E-02 


1.6E-01 


1.2 




Zl 


3.7E-02 


3.3E-02 


3.8E-02 


3.7E-02 


4.1E-02 


4.1E-02 


1.5 


Z2 


6.0E-02 


5.4E-02 


6.3E-02 


6.1E-02 


6.6E-02 


6.6E-02 


9.6E-02 



Note. — Absolute value of the noise cross-power-spectrum for the 145 
GHz detectors, averaged in the 0.1-2 Hz range. The numerical values have 
been normalized to the noise auto-power-spectrum of detector 145 Wl, which 
is 1.2 X IQ-^V^/Hz at the ADC input. 
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Fig. 21. — Flight-average of the auto power spectrum of the data time-stream, from channel 
145W1 (top line, red). The data are recorded at the ADC input. The total gain of the 
readout chain between the bolometer and the ADC input is 50000. The middle line (black) 
is the auto power spectrum of the time stream after subtraction of the signal estimated from 
the best fit map: this is an estimate of the instrumental noise. Comparing the two, it is 
evident that most of the CMB signal is encoded in the 0.05-1 Hz range. The lowest line 
(blue) is the cross correlation of the noise (time-stream - best fit map) signal from 145W1 
and 145W2. In the frequency range of interest for the CMB signal, the cross-correlated 
noise signal is at least one order of magnitude smaller than the auto-correlated signal. The 
correlated noise rises at very low frequencies, where atmospheric effects, thermal drifts, 
pendulations and scan-related microphonics can be present. At frequencies higher then 2 Hz 
the noise correlation remains well below the noise of the detector. 



where N is the noise covariance matrix, whose estimation is discussed in Sect. 6.5 above: 





N = (ntn*/) 



(28) 



This matrix becomes block diagonal assuming that there is no noise correlation between 
different detectors: 



Given the size of the problem, in the IT pipeline the map making normal equations (27) are 
solved by implementing a Fouricr-bascd preconditioned conjugate gradient iterative solver, 
which only needs to perform matrix to vector products. This approach critically depends 
on the assumption that the noise is stationary (Natoli et al. (2001)). In the NA pipeline a 
Jacobi solver is used. 

In practice, there are many further details. Since the solver works in Fourier space, 
the timeline must be continuous; bad or missing data chunks ("gaps") are replaced by a 
constrained realization of Gaussian noise (Hoffman & Ribak (1991)), designed to mimic the 
correct noise behavior at the gap's boundaries. The samples belonging to these chunks are 
then flagged and the solver is instructed not to take them into account, as they do not 
contain any useful data. The noise correlation functions (also called noise filters) are band 
pass filtered in the range 70 mHz - 14.5 Hz before being given to the map making code. 

Each detector has its own noise filter, which is kept constant over the full flight (i.e. 
data are reduced as if the noise were stationary within a single detector) for the IT pipeline, 
while the NA pipeline estimates noise separately for each (~ hour long) chunk of data. The 
NA approach treats the non stationarity of the data, at the cost of a reduced accuracy in 
the estimate of the noise filter. The differences in the resulting maps are neghgible. 

For the polarization maps, detectors forming PSB pairs are treated as independent 
detectors in the IT pipehne, while the difference of their signals is used by the NA pipehne 
(see Jones et al. (2005b) for details). 

The ROMA code employs about 20 minutes (and ~ 200 iterations) to produce a set of 
/, Q, U maps from an eight bolometer (full flight) timeline, when running a 128 processor 
job on an IBM SP machine featuring 450 MHz Power 3 processors. 



(nj4> = (n^'4> = {i^j)- 



(29) 
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7.2. Responsivity Calibration 

The first step in obtaining our final product maps of I, Q, U, is the relative calibration of 
all the detectors. As explained in §2.1 , we need a relative calibration accuracy of the order 
of ^ 2%. We start by producing single detector maps of I using one of the map-making 
procedures. We then proceed in two ways: in pixel space and in multipole space. 

7.2.1. Relative calibration from pixel-pixel correlation of different PSBs 

We carry out this analysis for the PSBs in the deep region. We use the 145 GHz Wi 
channel as the reference channel. The characteristics of these channels (beam, frequency 
response, noise, polarization efficiency) are so similar that their data can be directly com- 
pared, pixel to pixel. For channel j we scatter plot the li^j values (where i is the pixel index) 
versus the Ii;Wi values, and we fit the data to a straight line. We estimate the error on each 
as (7j/ Nj^i where aj is obtained by integrating the TOD noise power spectrum, and 
Nj^i is the number of observations of pixel i. The slope of the best fit line is the relative 
calibration factor TZj^Wi ■ Simulations show that this procedure results in an accurate (within 
1%) estimate of the relative calibration. The results for all channels are reported in table 8. 

We also study the effect of a mis-estimate of the noise in the responsivity estimate. This 
depends strongly on the signal to noise per pixel of the data. For 7' pixels, the signal to noise 
per pixel for CMB anisotropics is ~ 3. In these conditions, simple simulations show that 
in order to bias the relative calibration by more than 2% the estimate of the noise must be 
wrong by more than 20%. We conclude that this method provides a robust estimate of the 
relative calibration. The errors on the relative calibration results reported in table 8 include 
the effect of a conservative 20% misestimate of the errors per pixel. 

If, instead, the S/N per pixel is ^1, then a 20% misestimate of the noise induces 
^ 20% errors in the calibration. This happens when an absolute calibration is attempted 
by fitting lij versus the h^wMAP- Since the angular resolution of WMAP is worse than the 
one of BOS , and the noise per pixel is higher, the slope of the best fit line is significantly 
biased. Also, working in pixel space it is not trivial to take into account the difference in 
beam size and shape between BOS and WMAP . 
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Fig. 22. — Raw-data time-stream from channel 145W1. It is evident the scan-synchronous 
signal from the CMB dipole. The (red) line is the CMB dipole, measaured by WMAP and 
COBE, along this scan. There is clear agreement between the predicted CMB dipole and 
the observed scan synchronous signal. 



channel 


^ pixel 


^ spectrum 


1451^2 


0.956±0.028 


0.950 ±0.014 


145X1 


0.985±0.029 


0.980 ±0.016 


145X2 


0.749±0.022 


0.731 ±0.013 


145ri 


0.920±0.030 


0.932 ±0.020 


145F2 


0.912±0.029 


0.906 ±0.019 


145Z1 


0.674±0.021 


0.665 ±0.013 


145Z2 


0.407±0.023 


0.422 ±0.011 



Table 8: Relative calibration IZ of the PSB channels obtained from the pixel-pixel scatter 
plots with NSIDE=256 (second column, see §7.2.1) and from the cross-spectrum (third col- 
umn, see §7.2.2), using Wi as the reference channel. For 7?.^^^*™™, 1-a errors are used. For 
T^^*Wi-: also include in the error the possible bias due to a conservative ±20% error in the 
estimate of the noise. 



-60- 



1.2.2. Absolute and Relative calibration from ratios of cross-spectra 

We can compare the signals detected by different detectors observing the same sky in 
multipole space rather than in pixel space. A cross-power spectrum is very useful in this case, 
due to its unbiased nature: the noise properties of the two signals do not affect the results. 
Moreover, the use of the cross-power spectrum allows us to properly take into account in a 
simple way the effect of the different angular resolution, and map-making transfer functions. 
All these effects can be included in a function F^, which can be different for each channel k. 
Since we do see CMB fluctuations in all the BOS maps, this method can be used to calibrate 
the BOS detectors at 145, 245 and 345 GHz, and to find the absolute calibration factor by 
comparing BOS maps to WMAP maps. 

We expand the uncalibrated maps of all channels in spherical harmonics: 



14 = E 



i.m 



(30) 



We compute the cross power spectrum (a^^-F/ x a\^Fl) and compare it to the cross 
power spectrum {a^^^Fl x o^F/^^): the ratio is the relative calibration of channel k against 
the reference channel Wi. 



This is iterated for all j ^ k for consistency checks; the uncertainty on the result is 
estimated with a bootstrap method (see Polenta (2004) for details). The relative calibrations 
obtained in this way are reported in table 8. These numbers are consistent with the ones 
obtained in §7.2.1. We have also verified that using either of the two does not change the 
combined map. 

Having measured the relative calibrations for all channels of the same band, we can 
combine them into an uncalibrated optimal map of the sky, using the optimal map-making 
described in §7.1. 

We then estimate the absolute calibration of this map by comparing it to the WMAP map 
in the same region. The WMAP all-sky maps are calibrated with remarkable accuracy (0.5%, 
see Bennett et al. (2003)). So we can compare the angular power spectra measured by 
BOS and WMAP in the same sky region in order to measure the absolute responsivity S of 
BOS . 
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Assuming that the uncahbrated signal in BOS (in Volt) is related to the CMB tempera- 
ture measured by WMAP as IVmap = ^b/»5, we can estimate S from different combination 
of B03 and WMAP power spectra. We use the following one: 



(32) 



For the 145 GHz PSBs, we compute the cross-power spectrum of {dfmFfxa^Jf^^F^^^^) 
using optimal combined maps of all channels for both the experiments, and the cross-power 
spectrum of {af^^ F}^^ x af^^Fj^) using optimal maps obtained combining the W and 
X PSBs for the first map, and the Y and Z PSBs for the second one. Using this denominator 
in place of {af^^Ff x af^^^Ff) we avoid the need to remove a possible bias due to correlated 
noise. This choice, based on half of the detectors, does not degrade the precision of the 
result, which is already limited by cosmic variance. 

For the 145 GHz T map, we show in Figure 23 the measurement of S obtained from eq. 
32, i.e. the absolute calibration factor as a function of i. Its flatness confirms that we are 
properly correcting for beam, pixelization and noise effects. 

We confirm the robustness of the result by computing cahbration factors from the shal- 
low and deep surveys separately and by applying different weighting schemes. The final 
result for the 145 GHz map is 

cSi45 = (11.84 ± 0.20)V/K (33) 

where the 1-a error already incorporates the effect of the calibration uncertainty of WMAP . 
The absolute responsivity of a single channel k can be obtained from the values listed in 
table 8 as 5fc = STZk^Wi ■ 

The calibrations of 245 and 345 GHz maps are found in a similar way. We obtain 

^245 = (4.27 ± 0.20)V/K (34) 

and 

^345 = (3.04 ± Q.27)V/K (35) 



7.3. Differences in the two pipelines 



Before discussing the results of the analysis, we summarize in table 9 the most important 
differences between the two pipelines. Despite the many different choices, the final results 
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Fig. 23. — Absolute calibration factors Se based on WMAP measurement. The cross-power 
spectra are computed using the shallow integration region with uniform weight. The final 
value for S is the weighted average of the S^. 
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of the analysis are very consistent for both the maps (see §7.4 and following) and the power 
spectra (Jones et al. 2005; Piacentini et al. 2005; Montroy et al. 2005). 

7.4. Observations of the Galactic Plane 

l.Jf..!. Brightness Maps of the Galactic Plane 

Brightness maps of a section of the Galactic plane at 145, 245 and 345 GHz arc shown 
in Fig. 24, in comparison with other relevant maps of emission from dust and associated 
molecular gas. Several known Hll regions are evident in the maps. Particulary luminous 
are RCW38 (at RA=134.76°, dec=-47.47°) and 1RAS08576 (at RA=134.92°, dec=-43.75°). 
Diffuse emission is also evident. It is instructive to compare the B03 maps to monitors of 
different components of the ISM, like the IRAS/DIRBE maps for the dust continuum and 
the 115 GHz CO map of Dame et al. (2001) for the molecular gas (see bottom row of Fig.24). 
There is a good correlation between the BOS maps and the IRAS/DIRBE 3000 GHz map 
for most of the diffuse structure observed. A striking exception is the cloud evident in all 
B03 channels at RA=135.24° and dec=-44.80°: this is very dim in the IRAS map, while it 
is bright in the CO map and is also visible in the WMAP 94 GHz map: the signature of a 
cold dust cloud associated with molecular gas. 

The measured spectral flux density (SFD) of the main Galactic sources measured by 
BOS is reported in Fig. 25. The fluxes have been obtained integrating the brightness maps on 
disks centered on the sources. The diameter of the disks is chosen as the maximum between 
the apparent size of the source and twice the FWHM of the experiment beam. The error in 
the determination of the flux is dominated by the presence of diflusc emission, which must 
be subtracted to estimate the net flux of the source: this is less critical for RCW38, which 
is a relatively isolated source. This results in error bars larger than those from calibration 
and pointing jitter. 

In order to gain insight on the physical processes of emission operating in the sources, 
we combined the BOS FSD data with data from WMAP, DIRBE and IRAS (see Fig. 25). A 
combination of a power law dominating at low frequencies (~ i/^), and a dust-spectrum at 
high frequencies (~ v'^B{v,Td)) produces a good fit to the data. 

We find dust temperatures of (22.4 ± 0.9)i^, (19.6 ± 0.7)i^, and (13.3 ± 0.6)X for 
RCW38, IRAS08576, and the CO cloud respectively. The spectral index of the power law 
is (-0.5 ± 0.1), (-0.5 ± 0.3), and (-1.7 ± 1.3) for RCW38, IRAS08576, and the CO cloud 
respectively. 
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Fig. 24. — B03 survey of the Galactic plane. In the top row we display (left to right): 
B03 maps at 145, 245 and 345 GHz. The brightness units are MJy/sr. The conversion 
factors from MJy/sr to uiKcmb are given in table 3. In the bottom row we report, for 
comparison, maps of the same region from (left to right): the CO survey of Dame et al. 
(2001), the WMAP 94 GHz channel, and the IRAS 3000 GHz survey. All the maps have 
been convolved with a 6' FWHM gaussian filter. The equivalent FWHM resolution of the 
maps is indicated by the yellow circle on the bottom right of each panel. The hexagons 
locate RCW38 {dec 47.5°) and IRAS08576 {dec 43.8°). 
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Fig. 25. — Measured integrated fluxes for selected Galactic sources observed by B03 (dia- 
monds). The other data points are from WMAP (stars), from DIRBE (triangle) and from 
the IRAS/DIRBE map (square). The continuous line is the best fit SFD obtained as the 
sum of a power law at low frequencies (~ u"-) plus a thermal dust spectrum dominating at 
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A word of caution is necessary for the power-lax indices derived here. We are using 
data from experiments with different angular resolution. In particular the low frequency 
data have poorer angular resolution, so the measured flux can be contaminated by nearby 
sources entering these wider beams. For this reason the spectra shown in Fig. 25 could be 
increasingly contaminated (biased high) at low frequencies. This issue will be analyzed in a 
future publication. 

7.4.2. Polarization of RCW38 

RCW38 is the brightest source we observed. At 145 GHz the polarization of RCW38 is 
very low. We have carried out an analysis on the 145 GHz W, X, and Y PSBs, analyzing 53 
scans over the source, and assuming that there is a constant polarization over the size of the 
source. For each scan and for each of the 6 bolometers we fit a maximum amplitude of the 
detected signal. We then solve for Q and U of each scan using the detector signal differences 
W1-W2, X1-X2, Y1-Y2, taking into account the relative cahbration of the bolometers, the 
color corrections, and the orientation angles. The resulting average values for Q and U 
are Qrcw^s = (-0.046 ± 0.076)M Jy/sr and Urcwss = (0.054 ± 0.082)MJy/sr. Since the 
average brightness for RCW38 on the same scans is 3.7 MJy/sr, the upper limit for the 
polarization degree is Prcw38 < 6% (2-cr U.L.). 

7.5. Brightness Maps at High Galactic Latitude 

7.5.1. Maps of CMB Temperature Anisotropy at High Galactic Latitude 

The main product of this experiment is the 145 GHz map, which is shown as a large 
image in Fig. 26. The structure visible in the map with high S/N is CMB anisotropy. ISD 
is negligible in comparison to the CMB, as we show below. 

Maps of the shallow and deep regions at 145, 245 and 345 GHz are compared in Figs 
27 and 28 respectively. For a meaningful comparison, all maps have been filtered in the 
same way. Aggressive high-pass filtering of the time ordered data (cut-on at 7.5 times the 
scan frequency) was required for the 245 GHz and 345 GHz bolometers, in order to avoid 
artifacts due to scan synchronous noise. The brightness at 245 and 345 is due to a mixture 
of CMB and emission from interstellar dust (ISD), as evident from the comparison of the 
CMB-subtracted maps to the IRAS 3000 GHz map in the botton rows of Fig. 27. 

The temperature anisotropy of the CMB is measured with high S/N in the deep 145 
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Fig. 26. — 145 GHz I map from all the PSB bolometers of B03 . I is encoded in the 
false color scale in thermodynamic temperature units for a 2.725K blackbody. The pixel 
size is 3.4 arcmin. The data of this map will be made publicly available (together with 
a set of realistic simulations needed for quantitative analysis) at the BOS web servers: 
http://oberon.romal.infn.it/boomerang/b2k and http://cmb.phys.case.edu/boomerang . 
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Fig. 27. — The shallow and deep surveys of B03 . Top row: I maps from B03 . From left to 
right: 145 GHz; 245 GHz; 345 GHz. I is encoded in the false color scale in thermodynamic 
temperature units for a 2.725K blackbody. The pixel size is 3.4 arcmin. Bottom row, from 
left to right: IRAS/DIRBE image of the same region at 3000 GHz; B03 difference map 
obtained by subtracting the 145 GHz map from the 245 GHz map, in order to remove CMB 
anisotropy from the resulting map; BOS difference map obtained as 345 GHz map minus 145 
GHz map. The central region, where the noise is evidently lower, is the deep survey. In 
Fig. 28 we zoom on such region. The circles label the positions of known AGNs. To remove 
large-scale gradients, a high-pass filter at 7.5 times the scan frequency was applied on the 
time-ordered data: this filter is needed to remove artifacts due to scan synchronous noise in 
the 245 and 345 GHz maps. A gaussian filter with 7' FWHM was used as anti-aliasing. 
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Fig. 28. — The deep survey of BOS . Top row: I maps from BOS . From left to right: 145 
GHz; 245 GHz; 345 GHz. I is encoded in the false color scale in thermodynamic temperature 
units for a 2.725K blackbody. The pixel size is 3.4 arcmin. Bottom row, from left to 
right: IRAS/DIRBE image of the same region at 3000 GHz; BOS difference map obtained by 
subtracting the 145 GHz map from the 245 GHz map, in order to remove CMB anisotropy 
from the resulting map; BOS difference map obtained as 345 GHz map minus 145 GHz map. 
All the maps have been filtered in the same way. To remove large-scale gradients, a high- 
pass filter at 7.5 times the scan frequency was applied on the time-ordered data: this filter 
is needed to remove artifacts due to scan synchronous noise in the 245 and 345 GHz maps. 
A gaussian filter with 7' FWHM was used as anti-aliasing. 
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GHz map. Fig. 29 compares the 145 GHz deep map with maps of the same region obtained 
from BOOMERANG -98 at 150 GHz, and from WMAP at at 94 GHz. There is excellent 
morphological agreement between all three maps. 

To estimate more precisely the S/N of the 145 GHz deep map, we jackknife the data to 
provide two independent maps: one using data only from the first half of the deep survey 
(Dl) and one using data only from the second half of the deep survey (D2). We then form the 
sum and difference maps, (Dl+D2)/2 and (Dl-D2)/2. The sum map contains both signal 
and noise, the difference map only noise. For all maps we have used HEALPIX NS1DE=1024, 
i.e. a pixel side of 3.5'. Histograms of both maps, shown in Fig. 30, are gaussian-distributed. 
The rms in the sum and difference maps is cr+ = (93.2 ± 0.5)/ii^ and O- = (23.4 ± 0.2)iJ,K, 
respectively. The rms of the sky signal is thus V93.22 _ 23.42 ^ (99.2 ± 2.3)(xK, where 
the error now includes the uncertainty in absolute calibration. Realistic simulations of our 
observations, including all the details of scan speed, coverage, detectors noise, filtering, 
pixelization, etc. can be used to estimate the expected rms of CMB anisotropy in these 
observations, for the standard " concordance" model best fitting WMAP. We used the same 
simulations used for our spectral analysis (Jones et al. 2005; Piacentini et al. 2005; Montroy 
et al. 2005). The result is (94 ± 13)fiK (including cosmic variance). A similar result is 
obtained integrating the power spectrum of the " concordance" model over the B03 window 
function shown in Fig. 15. 

The distribution of the measured brightness in the deep 145 GHz map is accurately 
gaussian. The simplest gaussianity test is the evaluation of the skewness S3 and kurtosis S4 
of the pixel temperature distribution, computed as 



In Fig. 31 we compare the measured skewness {S3 = —0.062) and kurtosis {S4 = —0.053) to 
the distribution of S3 and S4 evaluated from our realistic simulations set. The agreement is 
very good: no hints for any deviation from gaussianity. This also confirms that systematic 
effects and foreground contamination are negligible in this map. 

The CMB structures mapped at 145 GHz are evident also in the 245 and 345 GHz maps 
(Fig. 28), even if with lower S/N, due to the higher background on the bolometers, to the 
lower number of detectors in these bands, and to the presence of increasing contamination 
from interstellar dust emission. 




(36) 
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item 


IT pipeline 


NA pipeline 


cross-polar response: 
time response: 
noise: 
Q,U map-making: 
map-making solver: 


full beam (table 4, col.4) 
in-flight, from cosmic rays (§6.1) 

estimated on the full flight 
all detectors treated individually 
Conjugate Gradients 


axial (table 4, col. 5) 
from lab measurements (§3.1) 
estimated in ~ hour long chunks 
PSB pairs differences 
Jacobi 



Table 9: Summary of the different choices made for the two data analysis pipehnes. 



WMAP W-bnnd: Stokes I B9B ISOGHz USGHz: Stokea I 

■J"" 300 uK_c 300 uK ''00 uK_t 

-300-20(^100 100 200 300 -300-200-100 100 200 300 -30(.l200- 100 100 2t:: :0( 




00 ' " 00 

HA [Des] KA [Des] 



Fig. 29. — Comparison of maps of the deep survey region from WMAP 94 GHz (left), the 
BOOMERANG 98 data at 145 GHz (center) and the BOS 145 GHz data (right). AH maps 
have 7' pixels. The improvement in the S/N of BOS is evident. A remarkable morphology 
agreement of the structures detected in three independent experiments is also evident. In all 
maps I is encoded in the false colors scale in thermodynamic temperature units for a 2.725K 
blackbody. 
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Fig. 30. — Comparison of the 145 GHz CMB I maps obtained from the two halves of the deep 
survey. The pixehzation of the map is 3.5 arcmin (HEALPIX NSIDE=1024). We histogram 
pixel temperatures resulting from the difference and sum maps. Only noise and systematic 
effects contribute to the former, while both signal and noise contribute to the latter. The 
two curves are labelled with their standard deviations. The standard deviation of the signal 
is thus V93.22 _ 23.42 = 90.2/i/s:. 
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Fig. 31. — The vertical lines represent the measured skewness (left panel) and kurtosis 
(right panel) of the pixel temperature distribution in the 145 GHz deep survey map. The 
histograms derive from realistic simulations of the measurement, assuming an underlying 
(gaussian) "concordance" model. 
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1.5.2. Foreground Contamination by Diffuse ISD Emission at High Latitude 

Our multiband maps allow us to carry out basic tests on the level of foreground con- 
tamination. We consider first the deep survey (carried out at Galactic latitudes \h\ > 25°). 
We use the SFD IRAS map D^ j at 100 /xm, corrected using DIRBE data (Schlegel et al. 
1999), as a template for ISD emission. We use the WMAP 94 GHz map Wij as a template 
for GMB anisotropy. Both maps have been sampled along the scans of the BOS , and then 
high-pass and low-pass filtered using the transfer function of the 145 GHz BOS detectors, in 
order to create the corresponding synthesized time-streams. These have been processed in 
the same way as the BOS detectors (see §7.5.1), to obtain the maps shown in Figs. 27, 28 
and 29. For this particular analysis we use a 14' pixelization, in order to avoid problems 
with the different beam size of the instruments, and carry out a linear fit in pixel space: 



where lij is one of the three BOS maps (145, 245, 345 GHz). The results of the fits are 
reported in table 10. The errors have been estimated by dividing the deep survey region 
into nine sub-regions, performing the fit for each subregion, and computing the standard 
deviation of the best fit parameters. In this way we account for variations in the properties 
of the foreground over the survey region. The A values are not reported since they are biased 
low, due to the relatively high noise of the WMAP map. Simulations show that, with the 
noise levels present here, this problem does not bias the estimates of the B values by more 
than a few %. This error is added in quadrature to the statistical error, which dominates 
the total error. 

The fact that the correlation coefficient for the WMAP template is always significant 
means that we have detected GMB anisotropy in all of our channels. The correlation coef- 
ficient for the IRAS template, instead, is statistically significant for the 245 and for the 345 
GHz survey. This means that the best fit values for B at 145 GHz provide only an upper 
limit for ISD contamination. As an example of a region where ISD fiuctuations are larger, we 
can consider the part of the Shallow Survey closer to the Galactic Plane: 90" < RA < 105°; 
—50° < dec < —40°. Proceeding as before, in this region we obtain the results reported in 
the lower part of table 10, which confirm the results found for the deep survey. 

We have carried out an analysis similar to the one just described, using the difference 
maps M345 = 7345 — 7i45 and M245 = I245 — I145 as " GMB-subtracted" maps. Since our 145 
GHz map is dominated by GMB anisotropy and features very low noise, M345 and Af2i5 
monitor all non-CMB brightness, and have lower noise than what we can obtain subtracting 
the 94 GHz WMAP map as a GMB template. The fit of equation 



lij = AWij + BDij + C 



(37) 



Mi,,- = S Aj + C 



(38) 
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Frequency 


R{A) 


R{B) 


B[fiK/{MJy/sr)] 




145 GHz 


0.634 


0.074 


(18 ±23) 


< 10 


245 GHz 


0.452 


0.173 


(49 ±21) 


12 ±5 


345 GHz 


0.331 


0.452 


(202 ± 18) 


50 ±5 


145 GHz 


0.588 


-0.085 


(-5 ±9) 


< 5 


245 GHz 


0.273 


0.143 


(37± 11) 


21 ±6 


345 GHz 


0.123 


0.380 


(183 ± 24) 


104 ± 14 



Table 10: Pixel-space fits of the B03 maps with a combination of a CMB template (from WMAP ) 
and a dust template (from IRAS), see eq. 37. For the three upper rows, the fit has been performed 
in the deep survey region ( 74" < RA < 90°; -50° < dec < -40° , - 32°, for a total of 2172 15' 
pixels). In this region, the rms fluctuation of the IRAS/DIRBE map is 0.25 MJy/sr. For the three 
lower rows, the fit has been done in the part of the Shallow Survey closer to the Galactic Plane 
(90° <RA< 105°; -50° < dec < -40°, ~ 23°, 2027 pixels of 15'). Here the rms fluctuation of 
the IRAS/DIRBE map is 0.57 MJy/sr. R{A) is the correlation coefficient for the CMB template, 
while R{B) is the correlation coefficient for the dust template. The last column gives the estimated 
brightness fluctuation due to the ISD component correlated to the IRAS/DIRBE map, in CMB 
temperature units. 
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gives results very similar to the fits to eq. 37 reported in table 10. 

The rms fluctuation of the dust template in the deep survey region is about 0.25 MJy/sr, 
so we get /:i^Tdust,rms = (12 ± 5)ijlKcmb at 245 GHz and /S.Tdust,vms = (50 ± 5)ijlKcmb at 345 
GHz. The upper hmit at 145 GHz is ATfiust,rms < ^OjiKcMB (1-c u.L). 

Finkbeiner et al. (1999) have studied the extrapolation to longer wavelengths of dust 
fluctuations detected by IRAS . Our results are consistent with their extrapolations. Their 
model 8 for the specific brightness of ISD is : 

Bdust{i^) oc qhu"^B{u,T,) + (1 - fiK'B{u,T,) (39) 

with ai = 1.67, Ti = 9.4X, ^2 = 2.7, T2 = 16.2K, fi = 0.0363, q = 13. The rms of the 
extrapolated maps in the deep survey region are 2, 11 and 47 fiKcMB in the 145, 245 and 
345 GHz bands respectively. So this model fits quite well our measurements (see also Fig. 
32). 

With no detection at 145 GHz, we use the spectrum of Model 8 to extrapolate - with 
all the necessary caveats - to 145 GHz. Integrating in our bands (specified in Fig. 9), and 
normalizing to our 345 GHz point, we get ATdust,rms = (2.9 ± 1.2)iiKcmb at 145 GHz. This 
value is much smaller than the CMB anisotropy at the same frequency (the total sky rms at 
145 GHz in our window function is ~ QOixK). 

The dust fiuctuations are non-gaussian. However, only very few structures deviate sig- 
nificantly from the rms stated above: these are the clouds visible in the 354-145 difference 
map at (RA,dec) = (88.5°, -48.3°) and (84.6°, -48.2°). Their size is < 1°, and their bright- 
ness, extrapolated at 145 GHz with Model 8, is ^ 20jjLK. 

Since the polarization of the diffuse cirrus at these frequencies is ^ 10% of the brightness 
(Pointhieu et al. 2005), we expect a polarized contribution from ISD ^ IfiK rms, which is 
small with respect to the expected CMB polarization [y/ {EE) ~ ?,^K rms in our beam and 
with 14' pixelization we are using for the analysis in this paragraph). 

The contamination estimates above depend on the assumption that ISD is well moni- 
tored by the IRAS/DIRBE template. They do not take into account the possible existence 
of a colder dust component, undetected by IRAS, and with a different angular distribution. 

Studying the residuals from the best fit of equations 37 and 38, we find that the 145 GHz 
residuals are consistent with instrumental noise, while there are, in the 245 and 345 GHz 
residuals, structures which could be either instrumental artifacts or real sky fluctuations, or 
a combination of the two. The SFD extrapolation to our frequencies fails to produce these 
features (the SFD dust temperature is in fact fairly constant, at 17-18 K, throughout the 
deep survey region). These structures are quite dim: the rms is ~ ^OjiK at 345 GHz and 
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~ lOjiK at 245 GHz. These numbers include the effect of instrumental noise and should 
be considered conservative upper limits for the sky fluctuations. Their extrapolation to 145 



Our deep survey region is representative of a fairly large fraction of the high latitude 
sky. To see this, we select increasing fractions of the sky by requiring that dust brightness, 
monitored by the IRAS/DIRBE map at 100 /xm, is equal or lower than a give threshold. We 
find that the best ~ 40% of the sky has the same rms fluctuation of dust brightness as in 
our deep survey, and that the best 75 % of the sky has a brightness fluctuation ^ 3 times 
larger than the one in our deep survey. 



A few point soiuxes are evident in the I maps. These AGNs have been used for testing 
the pointing reconstruction procedures as explained in §5.1 . The effect of the full population 
of resolved and unresolved AGNs as a contaminant in CMB anisotropy measurements is an 
important topic of discussion, in view of the forecasted ultra-sensitive surveys of the CMB 
anisotropy and polarization. The SED of Blazars is almost flat {Siv) ~ with ^ a ^ 0.3 
between few GHz and 100 GHz). At our three frequencies, bracketing the frequency of 
maximum brightness of CMB fluctuations, the contamination of the AGNs is thus expected 
to be minimal. 

Giommi & Colafrancesco (2004) have searched AGN catalogs for all the sources in the 
region observed by BOOMERanG . Prom their list of 54 sources, we find 8 AGNs in the 
deep survey region. Once we exclude the three brighter ones (which are evident in the map) 
their equivalent CMB temperature in a 9.5' beam is < 400/xir at 145 GHz. Since the same 
region is covered by about 5000 independent beams, it is evident that the presence of mm 
AGN emission cannot contaminate the morphology of the I maps. The integrated flux {S) 
and its fluctuations ((AS*)^) produced by all resolved and unresolved AGNs with differential 
LogN-LogS distribution ^ are given by 



GHz, using a reasonable dust spectrum like Model 8, gives AT^ 



rms ~ 



lOl^KcMB at 145 GHz. 



7.5.3. Foreground Contamination by Compact Sources 




(40) 



(41) 
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Fig. 32. — Data Points: IRAS-correlated dust fluctuations detected by BOS : the squares refer 
to tlie Deep Survey region (74° < RA < 90°; —50° < dec < —40°); tlie circles refer to the part 
of the Shallow Survey closer to the Galactic Plane (90° < RA < 105°; -50° < dec < -40°, 
2027 pixels of 15'). The continuous line is Model 8 of Finkbeiner et al. (1999), normalized 
at 345 GHz for the Deep Region. The dashed line is the same model, scaled to the shallow 
region via the ratio between the IRAS/DIRBE rms in the two regions. 
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When detected with an instrument with effective sohd angle these flux fluctuations pro- 
duce a rms signal equivalent to a CMB anisotropy ATagn given by 

{{AT^GNr) = T'aMB , j^^^^'^ ,,2 (42) 

Using the ^ and the SED of Giommi & Colafrancesco (2004), we get {{ASf) ~ 30, 25, 21 
Jy^/sr at 145, 245, 345 GHz respectively. Using the B03 beams to compute fl, we obtain 
^((ATagjv)^) = 5.8, 4.2, 6.0 fiK rms at 145, 245, 345 GHz respectively This is negligible 
with respect to the measured anisotropy of the map. 

AGN emission is usually polarized at a level p ^ 5%. Assuming random polarization 
directions and neglecting the natural dispersion in the p values, we get 

{{AQAGNy) = ^{{ATAGNf) (43) 

So we find ^/{{AQ^^^ = 0.20, 0.15, 0.21 fiK rms at 145, 245, 345 GHz respectively 

This is negligible with respect to the amplitude of the E-modc CMB polarization that 
wc expect in our maps (~ 3fiK rms for the 14' pixelization, and ~ AfiK rms for the 3.5' 
pixelization) , and to our statistical noise. Future experiments that seek to detect the much 
smaller B-mode polarization signal - particularly the lensing signal at small angular scales - 
will need to carefully account for the contribution to the B-mode signal by compact sources. 



7.6. Polarization of deep, high latitude maps at 145 GHz 

In Fig. 33 we present maps of the Stokes parameters Q and U at 145 GHz in the deep 
survey region. 

Unlike the 1 maps, here the CMB signal in each pixel is smaller than the noise. To 
estimate the signal to noise ratio (S/N) of the map, we have produced again two separate 
145 GHz maps, one from the first half (Dl) of the observations in the deep region, and another 
from the second half (D2). In Fig. 34 we compare the histogram of the Q (and U) data 
from the sum map (Dl+D2)/2 to the same histogram from the difference map (Dl-D2)/2. 
The difference histogram is very similar to the sum histogram, confirming that the signal 
is small with respect to the noise. For this reason we do not expect any CMB polarization 
structure to be visible "by eye" in the Q and U maps. The remarkable agreement of the sum 
and difference histograms in Fig. 34 indicates the stability of the system and the absence of 
systematic effects over the 6 days of measurement on the deep survey. 



-80- 



In Fig. 35 we plot the histogram of the measured polarization vector directions, com- 
puted as = arctan[?7j/Qj]. Since the region observed by B03 is large compared to the 
typical correlation scale (~ 1°) characteristic of CMB polarization, and since the two signals 
are dominated by uncorrelated noise, we naively expect a uniform distribution of the ctj. 
The data plotted in Fig. 35 approximately confirm this expectation. The sine-like deviation 
present in the data is due to the anisotropy of the noise. Even though the 8 polarized detec- 
tors cover all possible directions (in steps of 7r/8), the rotation of the focal plane with respect 
to the sky during the flight is small enough that for each direction of polarization only a 
few detectors have most of the statistical weight. If one of the bolometers is noisier than 
the others (see table 6), the excursions of the noise will be larger in that direction, and the 
measured polarization will be preferentially oriented in the same direction. This has been 
confirmed by repeating the same analysis on maps obtained from realistic noise simulations. 
The B03 data arc in full agreement with the distribution of the results of the simulations, 
which are plotted for comparison in Fig. 35. We can conclude that anisotropy effects are con- 
tained within ^ 10% and are well described by our simulations. These effects can thus be 
treated in a self consistent way by the map making and by the spectral estimation analysis 
and simulation procedures. 

From the sum and difference maps we compute the quantities CFq ^^y = (TQ,sum ~ ^Q,dif} 
and cr^^^y = o-ij^sum ~ '^u,diff- absence of systematics, these are estimators of the 

mean square polarized signal from the sky. From a gaussian fit to the Q histograms we 
find aq^sum = (42.5 ± 0.2)/xX, while crg^di// = (42.9 ± Q.2)^K. Analogously for U we get 
(^u,sum — (42.0 ± 0.2)//ir, while (Ju^diff — (41.8 ± 0.2)//K. These results do not allow one 
to extract the small sky signal (we expect ~ 4/ ^Krms for Q and U in the concordance 
model with 3.5' pixelization) , but are fully consistent with the level of noise of our receivers 
(see Fig. 36). 

In fact the uncertainty in (Jq^sky ^u,sky obtained from realistic simulations of the 
instrumental noise and of the observations is 25 ixK'^ (see Fig. 36). So we have (Jq sky = 
(-34 ± 2b)^iK\ and al 

sky ~ ^ 25)fiK'^. For the polarization degree we have — 
^Q,sky + ^u,sky — ^ 35)fiK'^. This sets an upper limit to the rms hnear polarization of 
the sky of Prms < T-^IJ'K (95% confidence level). This limit is close to the i/iRrms expected 
for the CMB polarization signal in the 3.5' pixelization. 

Power spectrum methods provide much more powerful tools to extract the sky polariza- 
tion from noisy maps, since it is much easier to separate different components according to 
their multipole content. In particular the effect of low- frequency noise and drift mostly af- 
fects the low multipoles £ ^ 50. While the rms analysis above does not separate the different 
multipoles, the spectral analysis does; artifacts at low multipoles can be removed, allowing a 
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much more sensitive probe for signal at higher multipoles. The apphcation of these methods 
to B03 is described in detail in Piacentini et al. (2005) and Montroy ct al. (2005), where it is 
shown that there is in fact a CMB polarization signal hidden in the Q and U maps described 
here. 



8. Conclusions 

The results of this flight of BOOMERanG demonstrate the effectiveness of balloon- 
borne 145 GHz PSBs in the measurement of CMB anisotropy and polarization, and in the 
survey of mm- wave emission from the interstellar medium. 

In our three frequency bands the long duration balloon platform offers an excellent trade- 
off between cost and performance. With bolometers cooled to 0.3K, observing through an 
uncoolcd telescope, the stratosphere at altitudes ^ 25 Km offers a very stable environment, 
and optimal loading conditions for these measurements. 

During the entire flight, the B03 cryogenics and electronics produced the correct environ- 
ment (in terms of operating temperature stability, radiative background level and stability, 
and electromagnetic disturbances) to operate the PSBs. These performed very closely to 
the theoretical limit, as described in §6.4 and 6.5. The attitude control system allowed us 
to reconstruct the pointing of the telescope with an accuracy of a few arcmin, sufficient for 
the main purpose of the experiment, i.e. the measurement of the power spectra of CMB 
anisotropy and polarization. 

BOS produced maps of the sky at 145, 245 and 345 GHz in two high-latitude regions 
and one region at low Galactic latitudes, in the Southern sky, producing a new survey of 
more than 1000 square degrees with ~ 10' resolution. At high latitudes degree and sub- 
degree scale anisotropy of the CMB is evident in the 145, 245 and 345 GHz maps. The rms 
fluctuation of the sky temperature measured from the 145 GHz map is perfectly consistent 
with the current ACDM model for the anisotropy of the CMB. The 145 GHz Stokes Q and 
U maps presented here allow a statistical (power spectrum) detection of CMB polarization 
at 145 GHz, as discussed in companion papers. 

The approach discussed here also shows the effectiveness of multi-band observations in 
monitoring the foreground signal. The Galactic and extragalactic foreground is negligible 

( ~ 10% compared to the cosmological signal both in anisotropy and in polarization) at 145 
GHz and at high Galactic latitudes. Based on the SFD maps, about 40 % of the sky has 
ISD contamination equal or lower than in the deep survey region studied here. Future B- 
modes searches, seeking a much smaller polarization signal, will be heavily affected by the 



-82- 



foregrounds measured here. 

The results presented here are comparable in sensitivity to those published using in- 
terferometric techniques at lower frequencies; however, this measurement is the first using 
bolometric detectors, which can be scaled to large numbers and high sensitivity for future 
investigations of CMB polarization. In the nearer term, the measurements presented here 
provide a glimpse at what the Planck HFI instrument will achieve at 145 GHz; using the 
same technology (PSB detectors), the HFI will reach over the whole sky roughly the same 
sensitivity per pixel as was achieved by BOS in the deep survey. 

We gratefully acknowledge support from the CIAR, CSA, and NSERC in Canada; Agen- 
zia Spaziale Italiana, University La Sapienza and Programma Nazionale Ricerche in Antar- 
tide in Italy; PPARC and the Leverhulme Trust in the UK; and NASA (awards NAG5-9251 
and NAG5-12723) and NSF (awards OPP-9980654 and OPP-0407592) in the USA. Addi- 
tional support for detector development was provided by GIT and JPL. GBN acknowledges 
support from a Sloan Foundation Fellowship, WCJ and TEM were partially supported by 
NASA GSRP Fellowships. Field, logistical, and flight support were supplied by USAP and 
NSBF; data recovery was particularly appreciated. This research used resources at NERSC, 
supported by the DOE under Contract No. DE-AC03-76SF00098, and the MacKenzie clus- 
ter at CITA, funded by the Canada Foundation for Innovation. We also thank the CASPUR 
(Rome-ITAUV^) computational facihties and the Apphed Cluster Computing Technologies 
Group at the Jet Propulsion Laboratory for computing time and technical support. Some 
of the results in this paper have been derived using the HEALPix package (Gorski, Hivon & 
Wandelt 1999). 

REFERENCES 

Abroe, M. E., Borrill, J., Ferreira, P.G., Hanany, S., Jaffe, A. H., Johnson, B. R., Lee, A. 
T., Rabii, B., Richards, P. L., Smoot, G. F., Stompor, R., Winant, C. D., Wu, J.H.P., 
2004, Ap.J., 605, 607 

Aliaga, A.M., Martinez-Gonzalez, E., Cayon, L., Argueso, F., Sanz, J.L., Barreiro, R.B., 

2003, New Astron. Rev., 47, 821 

Barkats, D., Bischoff, C, Farese, P., Fitzpatrick, L., Gaier, T., Gunderson, J., Hedman, M., 
Hyatt, L., McMahon, J., Samtleben, D., Staggs, S., Vanderlinde, K., & Winstein, B., 

2004, ApJ, submitted (astro-ph/0409380) 



-83- 



Battistelli, E.S., De Petris, M., Lamagna, L., Maoli, R., Melchiorri, F., Palladino, E., Savini, 
G., Mauskopf, P.D., Orlando, A., SPIE conference proceedings "Astronomical tele- 
scopes and instrumentation", astro-ph/0209180 

Bennett, C, Halpern, M., Hinshaw, G. , Jarosik, N., Kogut, A., Limonm M., Meyer, S., 
Page, L., Spergel, D., Tucker, G., WoUack, E., Wright, E., Barnes, C., Greason, M., 
Hill, R., Komatsu, E., Nolta, M., Odegard, N., Peirs, H., Verde, L., & Weiland, J. 
2003, ApJS, 148, 1 

Benoit, A., et al., 2003, A&A, 424, 571 

Bernardi, G., et al., 2003, Ap.J.Lett., 594, 5 

Bond, J. R., Jaffe, A. H., Knox, L., 1998, Phys.Rev., D57, 2117-2137 
Bond, J. R., Jaffe, A. H., Knox, L., 2000, Astrophys.J., 533, 19 

Boyle, L.A., Steinhardt, P.J., Turok, N., 2004, Phys.Rev. D70, 023504, hep-th/0307170 

Bridle, S.L. , Zehavi, I., Dckcl, A., Lahav, O., Hobson, M.P., Lasenby, A.N., 2001, 
Mon.Not.RoyAstron.Soc, 321, 333 

Brouw, W.N., Spoelstra, T.A.T., 1976, Astron. Ap. SuppL, 26, 129 

Bucher, M., Moodley, K., Turok, N., 2002, Phys.Rev., D66, 023528 

Caderni, N., et al., 1978, Phys.Rev.D., 17, 1901; Phys. Rev. D., 17, 1908 

Carretti, E., Tascone, R., Cortighoni, S., Monari, J., Orsini, M., 2001, New Astronomy, 6, 
173 

Carretti, E., McConnell, D., McClure-Grifiiths, N. M., Bernardi, G., Cortiglioni, S., Poppi, 
S., 2005, submitted to MNRAS, astro-ph/0503043 

Cartwright, J. K., Pearson, T. J., Readhead, A. C. S. , Shepherd, M. C., and Sievers, J. L., 
2005, Ap.J, in press, (astro-ph/0502174) 

Catalano, A., Conversi, L., De Gregori, S., De Petris, M., Lamagna, L., Maoli, R., Savini, G., 
BattisteUi, E.S., Orlando, A., 2004, submitted to New Astronomy, astro-ph/0405399 

Copeland, E.J., Kolb, E.W., Liddlc, A.R., and Lidsey, J.E. , 1993, Phys. Rev. Lett., 71, 219 
; 1993, Phys. Rev. D, 48, 2529. 



-84- 



S. Cortiglioni, G. Bernardi, E. Carretti, L. Casarini, S. Cecchini, C. Macculi, M. Ramponi, 
C. Sbarra, J. Monari, A. Orfei, M. Poloni, S. Poppi, G. Boella, S. Bonometto, L. 
Colombo, M. Gervasi, G. Sironi, M. Zannoni, M. Baralis, O.A. Peverini, R. Tascone, 
G. Virone, R. Fabbri, V. Natale, L. Nicastro, K-W.Ng, E.N. Vinyajkin, V.A. Razin, 
M.V. Sazhin, LA. Strukov, B. Negri, 2004, New Astronomy, 9, 297 

B. P. Grill, P. A. R. Ade, D. R. Artusa, R. S. Bhatia, J. J. Bock, A. Boscaleri, P. Gar- 
doni, S. E. Ghurch, K. Coble, P. de Bernardis, G. de Troia, P. Farese, K. M. 
Ganga, M. Giacometti, C. V. Haynes, E. Hivon, V. V. Hristov, A. lacoangeli, W. C. 
Jones, A. E. Lange, L. Martinis, S. Masi, P. V. Mason, P. D. Mauskopf, L. Miglio, 
T. Montroy, C. B. Ncttcrficld, C. G. Paine, E. Pascalc, F. Piaccntini, G. Polenta, 
F. Pongetti, G. Romeo, J. E. Ruhl, F. Scaramuzzi, D. Sforna, and A. D. Turner, 
2003, Ap.J.Supp.Ser., 148, 527 

T. M. Dame, D. Hartmann, and P. Thaddeus, 2001, Ap.J., 547, 792 

de Bernardis, P., De Gasperis, G., Masi, S., Vittorio, N., 1994, Ap.J. Letters, 433, LI 

P. de Bernardis, P. A. R. Ade, J. J. Bock, J. R. Bond, J. Borrill, A. Boscaleri, K. Coble, 

B. P. Grill, G. De Gasperis, P. C. Farese, P. G. Ferreira, K. Ganga, M. Giacometti, 
E. Hivon, V. V. Hristov, A. lacoangeli, A. H. Jaffe, A. E. Lange, L. Martinis, S. Masi, 
P. V. Mason, P. D. Mauskopf, A. Melchiorri, L. Miglio, T. Montroy, G. B. Netterfield, 

E. Pascale, F. Piacentini, D. Pogosyan, S. Prunet, S. Rao, G. Romeo, J. E. Ruhl, 

F. Scaramuzzi, D. Sforna, and N. Vittorio, 2000, Nature, 404, 955 

de Bernardis, P., Ade, P., Bock, J., Bond, J., Borrill, J., Boscaleri, A., Coble, K. , Contaldi, 

C, Grill, B., De Troia, G., Farese, P., Ganga, K., Giacometti, M. Hivon, E., Hristov, 
v., lacoangeli. A., Jaffe, A., Jones, W., Lange, A., Martinis, L. Masi, S., Mason, 
P., Mauskopf, P., Melchiorri, A. Montroy, T., Netterfield, G., Pascale, E., Piacentini, 

F. Pogosyan, D., Polenta, G., Pongetti, F., Prunet, S., Romeo, G., Ruhl, J., & 
Scaramuzzi, F. 2002, ApJ, 564, 559 

P. de Bernardis, P.A.R. Ade, A. Balbi, J.J. Bock, J.R. Bond, J. Borrill, A. Bosca- 
leri, P. Gabella, G.R. Contaldi, B.P. Grill, G. De Gasperis, A. De-Oliveira Costa, 

G. De Troia, K. Ganga, M. Giacometti, E. Hivon, V.V. Hristov, T. Kisner, A. la- 
coangeli, A.H. Jaffe, W.C. Jones, A.E. Lange, S. Masi, P. Mason, P.D. Mauskopf, C. 
MacTavish, A. Melchiorri, T. Montroy, F. Natl, P. Natoli, C.B. Netterfield, E. Pascale, 
F. Piacentini, D. Pogosyan, G. Polenta, S. Prunet, G. Romeo, S. Ricciardi, J.E. Ruhl, 
M. Tegmark, N. Vittorio, in astro-ph/0311396, lAU Symposium 216: Maps of the 
Cosmos. Sydney 14-17 July 2003 - ASP Conference Series 



-85- 



de Gasperis, G., Balbi, A., Cabella, P., Natoli, P., & Vittorio, N., 2005, A&A in press 
(astro-ph/0502142) 

De Troia, G., P. A. R. Ade, J. J. Bock, J. R. Bond, A. Boscaleri, C. R. Contaldi, B. P. Grill, 
P. de Bernardis, P. G. Ferreira, M. Giacometti, E. Hivon, V. V. Hristov, M. Kunz, 
A. E. Lange, S. Masi, P. D. Mauskopf, T. Montroy, P. Natoli, G. B. Netterfield, 
E. Pascale, F. Piacentini, G. Polenta, G. Romeo, and J. E. Ruhl, 2003, MNRAS, 343, 
284 

Delabrouille, J., Kaplan, J., & The Planck HFI Gonsortium. 2002, in AIP Gonf. Proc. 609, 
Astrophysical Polarized Backgrounds, ed. S. Gecchini, S. Gortiglioni, R. Sault, & G. 
Sbarra (New York: AIP), 135 

Djik, J., Van Dicpcnbcck, C. T. W., Maanders, E. J., Thurlings, L. F. G., 1974, I.E.E.E. 
Trans, on Antennas and Propagation, AP22, 513-520. 

Dodelson, S., Knox, L., 2000, Phys. Rev. Lett., 84, 3523 

Dore, O., Teyssier, R., Bouchet, F.R., Vibert, D. & Prunet, S., 2001, A&A, 374, 358; see 
also http : //ulysse . iap . f r/ cmbsof t/mapcumba/ 

Douspis, M., Bartlett, J. G., Blanchard, A., Le Dour, M. 2001, A&A, 368, 1 

Dragone, M., 1974, IEEE Trans. AP, 22, 472; 1982, IEEE Trans. AP, 30, 331. 

Efstathiou, G., Bond, J.R., 1999, MNRAS, 304, 75 

Farese, P., Dall'Oglio, G., Gundersen, J., Keating, B., Klawikowski, S., Knox, L., Levy, A., 
Lubin, P., O'Dell, G., Peel, A., Piccirillo, L., Ruhl, J. & Timbie, P. 2004, ApJ, 610, 
625 

Ferreira, P.G. & Jaffe, A.H., 2000, MNRAS, 312, 89 
Finkbeiner, D., Davis, M., & Schlegel, D., 1999, Ap.J., 524, 867 

Fischer, M. L. , Alsop, D. G. , Gheng, E. S. , Glapp, A. G. , Gottingham, D. A. , Gundersen, 
J. O. , Koch, T. C. , Kreysa, E. , Meinhold, P. R. , Lange, A. E. , Lubin, P. M. , 
Richards, P. L. and Smoot, G. F., 1992, Ap.J., 388, 242-252. 

Gordon, G., Lewis, A., 2003, Phys. Rev., D67, 123513 

Gervasi, M., Boella, G., Cavaliere, F., Grossetti, G., Passerini, A., Sironi, G., Tartari, A., 
Zannoni, M., astro- ph/030 1350 



-86- 



Giommi, P., & Colafrancesco, S., 2004, Astron. Astrophys., 414, 7. 

Gorski, K. M., Hivon, E., Wandelt, B, D., 1999, in Proceedings of the MPA/ESO Cosmology 
Conference "Evolution of Large-Scale Structure", eds. A.J. Banday, R.S. Sheth and 
L. Da Costa, PrintPartners Ipskamp, NL, pp. 37-42 (also astro-ph/9812350) 

Gundersen, J. O., et al.. New Astronomy Reviews, in press., astro-ph/0306010 

Guth, A., & Pi, S. Y., 1982, Phys. Rev. Lett., 49, 1110 

Hamilton, J.Ch, Benoit, A., the Archcops Collaboration, 2003, Comptes Rendus Physique 
Academie des Sciences, 4, 853, astro-ph/0310788 

Hanany, S., Matsumura, T., Johnson, B., Jones, T., Hull, J. R., Ma, K. B., 2003, IEEE 
Trans. Appl. Supercond., 13, 2128 

M. N. Hedman, D. Barkats, J. O. Gundersen, J. J. McMahon, S. T. Staggs, & B. Winstein, 
2002, ApJ, 573, L73 

Hoffman, Y. & Ribak, E. 1991, ApJ, 380, L5 

Holmes, W., Gildemeister, J. M. and Richards, P. L, 1998, Apphed Physics Letters, 72, 
2250. 

Holzapfel, W. L. , Wilbanks, T. M. , Ade, P. A. R. , Church, S. E. , Fischer, M. L. , Mauskopf, 
P. D. , Osgood , D. E. , Lange, A. E., 1997, Ap.J., 479, 17. 

Hu, W., White, M., 1997, New Astron., 2, 323 

Jaffe, A. H., Ade, P. A. , Balbi, A. , Bock, J. J. , Bond, J. R., Borrill, J., Boscaleri, A., Coble, 
K. B. P. Crill, P. de Bernardis, P. Farese, P. G. Ferreira, K. Ganga, M. Giacometti, 
S. Hanany, E. Hivon, V. V. Hristov, A. lacoangeli, A. E. Lange, A. T. Lee, L. Mar- 
tinis, S. Masi, P. D. Mauskopf, A. Melchiorri, T. Montroy, C. B. Netterfield, S. Oh, 
E. Pascalc, F. Piaccntini, D. Pogosyan, S. Prunet, B. Rabii, S. Rao, P. L. Richards, 
G. Romeo, J. E. Ruhl, F. Scaramuzzi, D. Sforna, G. F. Smoot, R. Stompor, C. D. 
Winant, and J. H. Wu, 2001, Phys. Rev. Letters, 86, 3475 

Johnson, B. R., Abroe, M. E., Ade, P., Bock, J., Borrill, J., Collins, J. S., Ferreira, P., 
Hanany, S., Jaffe, A. H., Jones, T., Lee, A. T., Levinson, L., Matsumura, T., Rabii, 

B. , Renbarger, T., Richards, P. L., Smoot, G. F., Stompor, R., Tran, H. T., & Winant, 

C. D. 2003, in The Cosmic Microwave Background and its Polarization, ed. S. Hanany 
& K.A. Olive (Amsterdam: Elsevier), in press (astro-ph/0308259) 



-87- 



Jones, R. J., 1953, J. Opt. Soc. Am, 43, 1. 

Jones, B., Bhatia, R.S., Bock, J.J., Lange, A.E., A Polarization Sensitive Bolometric Detec- 
tor for Observations of the Cosmic Microwave Background, 2003, SPIE, 4855, 227, 
astro-ph/0209132 

Jones, W. C. 2005, PhD thesis, Cahfornia Institute of Technology. 

Jones, B., et al., 2005, Ap.J., submitted. 

Jones, B., Montroy, T., Contaldi, C, Lange, A.E.L., Netterfield C. B., Ruhl, J.E., 2005, in 
preparation. 

Kaiser, N., 1983, MNRAS, 202, 1169 

R.E. Kalman, 1960, Journal of Basic Engineering, 82D, 35 

Kamionkowski, M. Kosowsky, A., & Stebbins, A. 1997, Phys. Rev. D, 55, 7368 

M. Kaplinghat, M. Chu, Z. Hainam, G.P. Holder, L. Knox, and C. Skordis, 2003, Ap.J., 583, 
24 

Keating B. G., O'Dell, C. W. , de Oliveira-Costa, A., Klawikowski, S., Stebor N., Piccirillo, 
L., Tegmark, M., Timbie, P. T., 2001, ApJ Letters, 560, LI 

Keating, B. G., O'Dell, C.W., Gundersen, J.O., Piccirillo, L., Stebor, N., Timbie, P.T., 2001, 
Ap.J.SuppL, 144, 1-20 

Keating, B. G., et al, http://www.astro.caltech.edu/' lgg/bicep/BICEP_color.pdf 

Kogut, A., Spergel, D. N., Barnes, C., Bennett, C. L., Halpern, M., Hinshaw, G., Jarosik, 
N., Limon, M. S., Meyer, S., Page. L., Tucker, G., WoUack, E., & Wright, E. L., 2003, 
ApJS, 148, 161 

Kolb, E.W., and Turner, M.S., The Early Universe, Addison- Welsey, 1990. 

Komatsu, E., Kogut, A., Nolta, M., Bennett, C. L., Halpern, M., Hinshaw, G., Jarosik, N., 
Limon, M., Meyer, S. S., Page, L., Spergel, D. N., Tucker, G. S., Verde, L., WoUack, 
E., Wright, E. L., 2003, Ap.J.SuppL, 148, 119 

Kovac, J., Leitch, E., Pryke, C., Carlstrom, J., Halverson, N., & Holzapfel, W. L. 2002, 
Nature, 420, 772 

Kuo, et al., 2004, Ap.J., 600, 32-51 



-88- 



Lange, A. E. , Preund, M. M. , Sato, S. , Hirao, T. , Matsumoto, T. , Watabe, T., 1994, 
Ap.J., 428, 384-392 

Lange, A., Ade, P.A.R., Bock, J., Bond, J., Borrill, J., Boscaleri, A., Coble, K. , Crill, 
B., de Bernardis, P., Parese, P., Perreira, P., Ganga, K., Giacometti, M., Hivon, E., 
Hristov, v., lacoangeli. A., Jaffe, A., Martinis, L., Masi, S., Mauskopf, P., Melchiorri, 
A. Montroy, T. Netterfield, C.B., Pascale, E., Piacentini, P., Pogosyan, D., Prunet, 
S., Rao, S., Romeo, G., Ruhl, J., Scaramuzzi, P., & Sforna, D., 2001, Phys.Rev., D 
63, 042001 

Leach, S.M., Liddle, A.R., 2003, Phys. Rev., D68, 123508 

Lee, C. , Ade, P. A. R. and Haynes, C. V., 1996, ESA SP-388, 81 

Leitch, E., et al, 2002, Nature, 420, 763. 

Leitch, E., Kovac, J., Halverson, N., Carlstrom, J., Pryke, C, & Smith, M. 2004, ApJ, 
submitted (astro-ph/0409357) 

Lewis, A., Bridle, S., 2002, Phys.Rev. D66, 103511 

Linde, A., 1983, Phys. Lett., B129, 177. 

Lubin, P., Smoot, G., 1981, Ap.J., 245, 1 

Mac Tavish C., et al., 2005, Ap.J., submitted. 

Masi, S., "Search for the Cosmic Background Polarization", in "Gamow Cosmology", R. 
Ruffini and F. Melchiorri eds, LXXXVI Course of the International School of Physics 
E. Permi, Varenna, 1984, 310-316, Soc. It. Pisica, Bologna. 

Masi, S., Aquilini, E., Cardoni, P., de Bernardis, P., Martinis, L., Raccanelli, A., Scaramuzzi, 
P., Sforna, D., 1998, Cryogenics, 38, 319 

Masi, S., Cardoni, P., de Bernardis, P., Piacentini, P., Raccanelli, A., Scaramuzzi, P., 1999, 
Cryogenics, 39, 217 

Masi, S., Ade, P. A. R., Bock, J. J., Boscaleri, A., de Bernardis, P., de Troia, G., di Stefano, 
G., Hristov, V. V., lacoangeli. A., Jones, W. C, Kisner, T., Lange, A. E., Mauskopf, 
P. D., Mac Tavish, C, Montroy, T., Netterfield, C. B., Pascale, E., Piacentini, P., 
Pongetti, P., Romeo, G., Ruhl, J. E., Torbet, E., & Watt, J. 2002, in AIP Conf. Proc. 
609, Astrophysical Polarized Backgrounds, ed. S. Cecchini, S. Cortiglioni, R. Sault, 
& C. Sbarra (New York: AIP), 122 



-89- 



Mather, J. C, 1982, Applied Optics, 21, 1125. 

Mauskopf, P.D., Bock, J. J., Del Castillo, H., Holzapfel, W.L., Lange, A.E., 1997, Applied 
Optics, 36, 4 

Montroy T., et al, 2005, Ap.J., submitted. 

Montroy T., 2003, Ph.D. Thesis, University of California at Santa Barbara 
( htt p : / / cmb . phy s . cwr u .edu/tom/thesis.ps) 

Mukhanov V.F., Chibisov G.V., 1981, JETP Lett., 33, 532-535 (1981, Zh. Eksp. Teor. Fiz., 
33, 549-553) see also astro-ph/0303077 

Nanos, CP., 1979, Ap.J., 232, 241. 

Natoh P., de Gasperis, C, Gheller, C. & Vittorio, N., 2001, A&A, 372, 346 

Natoh, P., Marinucci, D., Cabella, P., de Gasperis, G., & Vittorio, N. 2002, A&A, 383, 1100 

Netterfield, C. B., et al., 1997, Ap.J., 474, 47 

Netterfield, C, Ade, P., Bock, J., Bond, J., Borrill, J., Boscaleri, A., Coble, K. , Contaldi, C, 
Grill, B., de Bernardis, P., Farese, P., Ganga, K., Giacometti, M. Hivon, E., Hristov, 
v., lacoangeli. A., Jaffe, A., Jones, W., Lange, A., Martinis, L. Masi, S., Mason, P., 
Mauskopf, P., Melchiorri, A. Montroy, T., Pascale, E., Piacentini, F. Pogosyan, D. 
Pongetti, F., Prunet, S.,Romeo, G., Ruhl, J., & Scaramuzzi, F. 2002, ApJ, 571, 604 

Padin, S., et al., 2002, Pub. Astron. Soc. of Pacific, 114, 83. 

Partridge, B., et al, 1988, Nature, 331, 146. 

Peiris, H.V., et al., 2003, Ap.J.Suppl., 148, 213 

F. Piacentini, P. A. R. Ade, R. S. Bhatia, J. J. Bock, A. Boscaleri, P. Cardoni, B. P. Grill, 
P. de Bernardis, H. Del Gastillo, G. De Troia, P. Farese, M. Giacometti, E. F. Hivon, 
V. V. Hristov, A. lacoangeU, A. E. Lange, S. Masi, P. D. Mauskopf, L. Miglio, G. B. 
Netterfield, P. Palangio, E. Pascale, A. Raccanelli, S. Rao, G. Romeo, J. Ruhl, and 
F. Scaramuzzi, 2002, Ap.J.Suppl.Ser., 138, 315 

Piacentini F., et al., 2005, Ap.J., submitted. 

Piccirillo, L., Ade, P. A. R., Bock, J. J., Bowden, M., Ghurch, S. W., Ganga, K., Gear, 
W. K., Hinderks, J., Keating, B. G., Lange, A. E., Maffei, B., Mallie, O., Melhuish, 
S. J., Murphy, J. A., Pisano, G., Rusholme, B., Taylor, A., & Thompson, K. 2002, 



-90- 



in AIP Conf. Proc. 609, Astrophysical Polarized Backgrounds, ed. S. Cecchini, S. 
Cortiglioni, R. Sault, & C. Sbarra (New York: AIP), 159 

Pointhieu, N., et al., 2005, A&A submitted, astro-ph/0501427 

Polenta, G., Ade, P.A.R., Bock, J., Bond, J., Borrill, J., Boscaleri, A., Coble, K. , Contaldi, 
C, Grill, B., de Bernardis, P., Farese, P., Ganga, K., Giacometti, M. Hivon, E., 
Hristov, v., lacoangeli. A., Jaffe, A., Jones, W., Lange, A., Martinis, L. Masi, S., 
Mason, P., Mauskopf, P., Melchiorri, A. Montroy, T., Pascale, E., Piacentini, F. 
Pogosyan, D. Pongetti, F., Prunet, S., Romeo, G., Ruhl, J., & Scaramuzzi, F., 2002, 
Ap.J., 572, L27 

Polenta, G., 2004, Ph.D. Thesis, University of Rome La Sapienza. 

Prunet, S., et al., 2001, proc. of the MPA/ESO/MPA conference "Mining the Sky" [astro- 
ph/0101073] 

Richards P.L., 1994, J. Appl. Phys. 76, 1. 

Readhead, A. G. S., Myers, S. T., Pearson, T. J., Sievers, J. L., Mason, B. S., Gontaldi, G. R., 
Bond, J. R., Altamirano, P., Bustos, R., Achermann, G., Bronfman, L., Garlstrom, J., 
Gartwright, J. K., Gasassus, S., Dickinson, G., Kovac, J., Leitch, E., May, J., Padin, 
S., Pogosyan, D., Pospieszalski, M., Pryke, G., Reeves, R., Shepherd, M., & Torres, 
S. 2004, Science, 306, 836 (astro-ph/0409569) 

Rees, M., 1968, Ap.J.Lett., 53, LI 

G. Romeo, P. de Bernardis, G. di Stefano, S. Masi, F. Piacentini, F. Pongetti, and S. Rao 
2002, In AIP Conf. Proc. 616: Experimental Cosmology at Millimetre Wavelengths, 
59. 

Ruhl, J. E., Ade, P., Bock, J., Bond, J., Borrill, J., Boscaleri, A., Gontaldi, G., Grill, B., 
de Bernardis, P., De Troia, G., Ganga, K., Giacometti, M. Hivon, E., Hristov, V., 
lacoangeh. A., Jaffe, A., Jones, W., Lange, A., Masi, S., Mason, P., Mauskopf, P., 
Melchiorri, A., Montroy, T., Netterfield, G., Pascale, E., Piacentini, F. Pogosyan, D. 
, Polenta, G., Prunet, S., Romeo, G., 2003, ApJ, 599, 786 

Savage, R., ct al., 2004, Mon.Not.RoyAstron.Soc, 349, 973 

Schlegel D.J. et al. 1999, Ap.J. 500, 525. 

Seljak, U., Pen, U., Turok, N., 1997, Phys. Rev. Lett., 79, 1615 



-91 - 



Song, Y.S., Knox, L., 2003, Phys. Rev., D68, 043518 
Spergel, D.N., et al., 2003, Astrophys.J.SuppL, 148, 175 
Steinhardt, P.J., Turok, N., 2002, Science, 296, 1436 

Subrahmanyan, R., Kesteven, M. J., Ekers, R. D., Sinclair, M., & Silk, J. 2002, MNRAS, in 
press (astro-pli/0002467) 

Tegmark, M., Zaldarriaga, M. 2000a, ApJ, 544, 30 

Tegmark, M., Zaldarriaga, M., 2000, Phys. Rev. Lett., 85, 2240 

Tegmark, M., et al., 2004, Phys.Rev. D69, 103501 

Turner, M.S., 1993, Phys. Rev. Lett., 71, 3502; 1993, Phys. Rev. D., 48, 5539. 

Villa, F., Mandolesi, N., Bersanelli, M., Butler, R. C, Burigana, C, Mennella, A., Morgantc, 
G., Sandri, M., & the LFI Consortium. 2002, in AIP Conf. Proc. 609, Astrophysical 
Polarized Backgrounds, ed. S. Cecchini, S. Cortiglioni, R. Sault, & C. Sbarra (New- 
York: AIP), 144 

WoUack, et al., 1997, Ap.J., 476, 440 

Yamada, Y., Mitsuishi, A., and Yoshinaga, H., 1962, J. Opt. Soc. Am., 52, 17. 
Yun, M., et al., 2003, Proc. SPIE, 4855, 136. 
Zaldarriaga, M., astro-ph/0305272. 

Zaldarriaga, M., Seljak, U., 1998, Phys. Rev., D58, 023003 
Zaldarriaga, M., Seljak, U. 1997, Phys. Rev. D, 55, 1830 



This preprint was prepared with the A AS IM^jX macros v5.2. 



-92- 




Fig. 33. — Polarization in the deep survey. Left: Q map from the 145 GHz PSB bolometers; 
right: U map from the same bolometers. Note that the stretch of the false color scale is a 
third of the one in the T maps. These maps are dominated by instrumental noise. 
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Fig. 34. — Comparison of the 145 GHz maps of Stokes Q (top) and U (bottom) taken in the 
two halves of the deep survey. The pixehzation has NSIDE=1024 (3.5' pixels). We histogram 
pixel values resulting from the sum (heavy line) and difference (light line) maps. The sum 
and difference histograms are very similar: at the level visible here there is no indication of 
systematic effects that could arise from changes of observation conditions during the survey. 
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Fig. 35. — Histogram of the directions of the polarization vectors, computed as a = 
arctan[L'"/Q], in the deep survey at 145 GHz (blue line). The ~ 15% sinusoidal devia- 
tion from uniformity is due to unequalized noise of the different PSBs. This is confirmed by 
the analysis of maps obtained from realistic noise simulations (black data points with l-o" 
error bars). 
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Fig. 36. — Histogram of tlic quantity (JQ^sum ~ ^Q,diff simulations of the B2K observa- 
tions of the Deep Region. The simulations used the concordance model to generate CMB 
maps and realistic noise simulations for all the 145 GHz detectors. The measured values of 
'^Q,sum - ^Q,diff and (7^_^„^ - a^aiff are plotted as vertical dashed hues. 



